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Abstract 

by 

MARY ELIZABETH CONVERY 



MiniMax (Fermilab T-864) was a small test/experiment at the Tevatron de- 
signed to search for disoriented chiral condensates (DCC) in the forward 
direction. 

Relativistic quantum field theory treats the vacuum as a medium, with 
bulk properties characterized by long-range order parameters. This has led 
to suggestions that regions of "disoriented vacuum" might be formed in high- 
energy collision processes. In particular, the approximate chiral symmetry 
of QCD could lead to regions of vacuum which have chiral order parameters 
disoriented to directions which have non-zero isospin, i.e. disoriented chiral 
condensates. A signature of DCC is the resulting distribution of the fraction 
of produced pions which are neutral. 

The MiniMax detector at the CO collision region of the Tevatron was a 
telescope of 24 multi-wire proportional chambers (MWPC's) with a lead con- 
verter behind the eighth MWPC, allowing the detection of charged particles 
and photon conversions in an acceptance approximately a circle of radius 0.6 
in pseudorapidity-azimuthal-angle space, centered on pseudorapidity 77 ft; 4. 
An electromagnetic calorimeter was located behind the MWPC telescope, 
and hadronic calorimeters and scintillator were located in the upstream anti- 
proton direction to tag diffractive events. 

The use of standard Monte Carlo simulations for high-energy collisions 
of elementary particles (PYTHIA) and for interactions of particles in the 
detector (GEANT) is described, along with the simulation created by the 
MiniMax Collaboration to generate DCC domains. 

A description of the data analysis software is given, including detailed 
studies of its performance on data from the simulations. 
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A set of robust observables is derived. These are insensitive to many 
efficiencies and to the details of the modehng of the parent pion production 
mechanisms, yet have distinguishable values for DCC and generic charged- 
neutral distributions. Simulations show that the robust observables are in- 
sensitive to detector efficiencies and to systematic errors in the data analysis 
software. 

The resulting values for robust observables for approximately 1.5 x 10^ 
events are shown to be consistent with production by only generic mecha- 
nisms. Results from samples of diffractive-tagged events and of high-multiplicity 
events also show no evidence for DCC. 
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Chapter 1 
Introduction 



1.1 Physics motivation 

1.1.1 Disoriented Chiral Condensates 

The Lagrangian of QCD (Quantum Chromodynamics the theory of strong 
interactions) for two quarks (u and d) has isospin symmetry SU(2). In the 
hmit that the quarks are massless, the Lagrangian also has chiral symmetry 
SU(2)ixSU(2)j:j. That is, if we write the isospin doublet 




in terms of left- and right-handed fields, 



^ = -(l-75)^ + -(l + 75)^ = ^L + ^R, (1.1) 

where 75 is the usual product of gamma matrices of Dirac theory, then the 
Lagrangian is symmetric under \E'i \E'/j. The spontaneous breaking of 
this chiral symmetry in the QCD ground state is accompanied by a massless 
Goldstone boson 0, the pion. In the real world, the quarks are light but 
not massless, and the QCD Lagrangian has an approximate chiral symmetry 
which is explicitly broken, giving the pion a small mass. 

Spontaneous chiral symmetry breaking is often described by the linear 
sigma model [^, in which there is a scalar field a of isospin and a vector 
pion field vf of total isospin 1. The Lagrangian can be written as 

^ = kd.nd^^ + d.ad^'a) - ^(vf^ + - fjf, (1.2) 



1 



2 



and the potential is minimized for cr^ + vf^ = /^^. The symmetry is broken 
when a particular minimum is chosen. In order to break chiral symmetry, 
but not isospin symmetry, the minimum chosen is that for which the a field 
acquires a vacuum expectation value: (a) = /^r, (vr) = 0. A term can be 
added to the Lagrangian in order to break the chiral symmetry explicitly, 
such as fT^m^a or — |m^7r^, but it is unclear experimentally which of these 
terms is realized in nature 

The ordinary vacuum has chiral order parameter in the a direction (it 
has no isospin). A disoriented chiral condensate (DCC) is a piece of vacuum 
which is disoriented from the a direction to a direction with vf components. 
When the DCC domain makes contact with the ordinary vacuum, it coher- 
ently radiates pions with isospin determined by the direction of disorientation 
in order to restore the a direction. For example, if the disorientation were 
in the 7r° direction, vr'^'s would be emitted. There has been much theoretical 
work done recently on DCC and other mechanisms for coherent semiclassi- 
cal radiation of pions in high-energy collisions of hadrons and of heavy ions 

The distinctive signature of DCC is that the pions produced when the 
DCC domain makes contact with the outside vacuum have a neutral fraction 
/ = Nt,o/{N^+ + A^^- + Nt,o) distributed according to 

P{f)df = ^^^df (1.3) 

in the limit of large numbers of pions. This neutral fraction distribution 
is common to some other mechanisms which produce coherent final states 



r^-[|T4|, [|T^-[|T9]. The proof of this has been given in terms of quantum 



mechanical coherent state arguments , but can be seen easily with a geo- 



metrical argument. The vacuum condensate has equal probability of having 
an order parameter oriented in any direction in {a, tt) space. If we define 6 as 
the polar angle relative to the vr'^-axis, quantum-mechanical arguments give 
/ = cos^e, and Pif)df = P{cos^ 6) d{cose). We have df/d^cosO) = 2//, 
so that P{f)df = {l/2y/J) df . In generic particle production, that is, pro- 
duction by (observed) mechanisms other than DCC, producing a pion of any 
given charge is equally likely due to isospin symmetry, so that / is binomially 
distributed with mean (/) = !. Note that the l/{2^) distribut ion also has 
a mean of (/) = |, but that the distributions are quite different; in particu- 
lar, the probability of producing pions with a small or large neutral fraction 
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is much higher for DCC than for generic production (Fig. |1 . 1| ) . 



It is conjectured that such a condensate may be formed in hadron-hadron 
coUisions with high transverse energy and a large multiphcity (number of 
particles produced). As the collision debris expands outward at almost the 
speed of light, it may form a hot, thin shell, the cool interior of which would 
be separated from the outside vacuum and could conceivably have a chiral 
order parameter which is disoriented from the a direction. When the shell 
hadronizes, the interior condensate makes contact with the outside vacuum 
and radiates pions. 

1.1.2 Centauro/anti-Centauro 

Cosmic ray experiments have found evidence for hadronic events which can be 
interpreted as having an anomalously large or small fraction of pions which 
are neutral and may therefore be related to DCC. Centauro events have 
been observed in emulsion chambers by the Chacaltaya-Pamir Collaboration 



2C, 21| and are characterized by a large number (~ 100) of charged particles 
and almost no electromagnetic energy, which implies no vr^'s, since their 
immediate decay produces two photons. The center-of-mass (cm) energies 
are on the order of a TeV or larger. Further interpretation of the events is 
controversial. What is measured can be related to the transverse momentum, 
e.g., (pt) = 0.35±0.15GeV for Centauro I, but the gamma ray inelasticity 
is not known. The value is usually quoted as kj ~ 0.2 — 0.4 with the 
lower and upper values preferred for nucleons and pions, respectively. Most 
analyses assume that the hadrons produced in Centauros events are nucleons. 
However, if we want to interpret Centauros as being related to DCC, the 
hadrons should be pions, and therefore it is possible that the pt is low {{pr) ~ 
0.875 ± 0.375 GeV with a large systematic uncertainty), which is referred to 
as Chiron behavior. It has been suggested (Ref. based on data from Ref. 
[PT|]) that these events may be diffractive.Q Taking the view that Centauros 
are diffractive fireballs recoiling against a proton or anti-proton, but that 

-'^DifFractive processes are thought of as involving the exchange of a colorless object 
called a pomeron, but are not well understood. In single diffraction, either the incoming 
proton or anti-proton is dissociated, while the other remains intact and typically has a 
low transverse momentum and a longitudinal momentum almost that of the initial beam. 
Both hadrons are dissociated in double diffraction. 
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pions rather than nucleons are produced, boosting these events to the lab 
frame of Fermilab Tevatron colhsions, Centauros would be expected to occur 
at ~ 3.5 — 4.5.0 (As will be discussed in Sec. the acceptance of the 
MiniMax detector covers precisely this region.) 

Events with a large neutral fraction are referred to as anti-Centauros. 
Such cosmic ray events have been reported by the JACEE Collaboration 
[p^ , which uses balloon-borne detectors. An example is the event shown in 
Fig. |1.2| . The leading cluster contains about 32 photons and only 1 charged 
particle. A possibly distinct cluster has about 54 photons and 17 charged 
particles. The collision occurred within the detector with a cm energy greater 
than 200 GeV and Chiron behavior was exhibited. 

However, the interpretation of these cosmic ray events is controversial, 
and observations of this type of event under controlled conditions, such as 
in a collider environment, where they could be better understood, would 
be very useful. Several collider experiments have unsuccessfully looked for 
large Centauro-like domains. CDF has conducted such a search at the Teva- 
tron {^/s = 1.8 TeV), and has reported no evidence for Centauros or anti- 
Centauros . However, the CDF search only covered the region | | ^ 3 
and only looked at particles with pt ^ 400 MeV, and therefore might not 
be expected to observe such phenomena. The UA5 Collaboration ruled out 
large Centauro domains going out to larger t] EO, |2^, but at cm energies 



no greater than 900 GeV, whereas the necessary cm energy is expected to 
be larger. UAl also conducted a search at relatively low energies and in the 
central region, and found no evidence [2^ . 



1.1.3 Multiplicity distributions 

Very little is known about particle production in the forward direction be- 
cause this region has not been well studied. The distribution of charged 



particles in lego space has been reported by experiments such as CDF p9 | 



^We work in what is called "lego space" where is the azimuthal angle and rj is the 
pseudorapidity which is defined in terms of 9, the angle from the beam axis, as = 
— In tan I . Therefore, the region explored by central detectors which look transverse to 
the beam covers small 77, while tj gets infinitely large in the beam direction. An advantage 
to using the pseudorapidity is that 77 has a simple transformation under boosts in the 
beam direction. Also, the density of particles produced in a collision is uniform over a 
large region of lego space. 



5 



and UAl for the central region (| 77 1 < 3). Measurements of dNch/drj were 
made for larger t] by UA5 |2^, and by P238 at the CERN SppS collider 
for proton-anti-proton collisions at center-of-mass energy = 630 GeV. 
The distribution of photons dN^/drj in the forward direction is even less 

, H. MiniMax, 



well known; UA5 observed photons at ^/s < 900 GeV | 
which was able to observe both charged particles and photons in the re- 
gion 3.3 < 7] < 4.5 at y/s = 1.8 TeV, therefore had the potential to make 
measurements in a previously unexplored region. 

Figure |L3| shows a plot of dNch/dr], averaged over typical collisions at 
Tevatron energies, from the event generator PYTHIA. The distribution in 
the central region is taken from measured values, and these are extrapolated 
to the forward region. Note that the mean number of charged particles 
varies by only about 0.5 (about 13%) in the region \ \ < 4. Since there is 
no preferred azimuthal direction, the distribution in is also uniform. 



1.2 Conceptual design of MiniMax 

With the primary goal to search for DCC, the MiniMax detector should be 
able to observe both charged particles and photons simultaneously. If we 
take seriously the interpretation of Centauros and anti-Centauros given in 
Sec. p..l.2| , the detector should cover the forward region 3.5 < r] < 4.5, 
and be sensitive to low-pr particles. The smallest coverage which would be 
expected to be sufficient for observing such events is ArjAcj) = 1. Further 
considerations include the restricted area around the beampipe due to the 
main ring and the floor, and the lack of funding for the experiment. In or- 
der to achieve "minimal maximum acceptance" (hence the name MiniMax), 
the detector was designed as a telescope of multi-wire proportional chambers 
(MWPC's) along the beampipe, with converter inside the telescope. Thus, 
charged particles can be observed in the chambers before and after the con- 
verter, and photon conversions in the chambers behind the converter. 

In the absence of a magnetic fleld, the energy of the particles can only 
be determined using calorimetry. For this reason, and to observe photons 
which do not convert inside the telescope, an electromagnetic calorimeter 
was placed behind the MWPC telescope. 

If Centauros are diffractive and are related to DCC, then DCC might be 
more likely to be produced in diffractive interactions than in non-diffractive 
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ones. In order to test this conjecture, diffractive events must be identified. 
Scintillator and hadronic calorimetry were used to detect leading anti-protons 
and anti-neutrons from diffractive events. 

A picture of the original design for the MiniMax detector is shown in Fig. 

L|. 

A brief note on the similarity between MiniMax and UA5 is in order, 
although the details of the UA5 experiment did not play a role in the design 
of the MiniMax detector. UA5 was an experiment at the CERN proton- 
anti-proton collider (SppS) and ran from 1981-1982 and in 1985, during 
which time collisions at ^/s = 546, 200, and 900 GeV were studied. The 
detector consisted of two streamer chambers, 600 x 125 x 50 cm^, on op- 
posite sides of the beampipe, in which charged particles left tracks which 
were photographed for analysis by well developed techniques. For the runs 
at y/s = 546 GeV, photons were detected through conversions in lead-glass 
plates, approximately 1 Xq thick, located inside the streamer chambers near 
the sides. In later runs, a lead converter plate (2 mm of lead supported by 
an aluminum box with walls 1 mm thick), placed between the beam pipe and 
the upper streamer chamber, was used instead. The beampipe was elliptical 
with mean dimensions 6 x 15.2 cm^. Combining this transverse distance of 
the chambers from the beam axis with the length of the chambers, UA5 was 
able to observe particles at pseudorapidities up to \ri\ ^5. UA5 published 
results on many important studies, including multiplicity and pseudorapidity 
distributions and correlations for charged particles and photons, strangeness 
production (e.g. K's, A's, and S's, observed through decays), diffractive 
dissociation, and a search for Centauros. Reference |^ is a comprehensive 
report on these studies at ^/s = 546 GeV. The Monte Carlo simulations cre- 
ated and used by UA5, which include non-diffractive and diffractive event 
generators, along with a Centauro generator, are described in Ref . . 




Figure 1.1: Binomial and DCC neutral fraction distributions. 
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Figure 1.2: JACEE event 
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Figure 1.3: Charged multiplicity distribution dNch/drj for non-single diffrac- 
tive inelastic PYTHIA events at ^/s = 1.8 TeV. 
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Figure 1.4: Original design of the MiniMax detector. 



Chapter 2 
MiniMax 



2.1 The MiniMax environment 

The MiniMax detector was located at the CO colhsion region of the Fermilab 
Tevatron, which colhdes protons and anti-protons at cm energy ^/s = 1.8 



TeV. A sketch of the detector is shown in Fig. p.l| . The coordinate system was 
taken to be left-handed with positive z in the downstream proton direction 
and positive y upwards. The region included a pit which was 18 in deep, 60 in 
wide, and 170 in long, with the bottom centered at {x, y, z) = (0, —29 in, 1 in) 
relative to the nominal collision point defined as (0,0,0). (Note that the 
actual collision point as determined by the mean distance of closest approach 
between tracks is found to be at 2 ~ 7 in, and within 0.5 in of x = ?/ = 0. 



The pointing of the tracks will be discussed in more detail in Sec. ^]3|.) 

The Tevatron beampipe in the neighborhood of the CO collision point 
during early stages of the experiment (from late 1993 to early 1995) was a Gin- 
diameter, 0.0625 in-thick Al pipe, with a 2 in-thick Al flange at the transition 
to a 2.5 in-diameter, 0.035 in-thick Al pipe and an abort pipe of the same 
dimensions. The detailed dimensions and locations of the pipe segments are 



given in Table |2.1| . Secondaries from interactions of collision primaries in the 
flange dominated the particles observed in the forward region. 

A new beampipe was designed to minimize the number of interactions in 
the pipe which produced background hits in the detector. This was accom- 
plished by steps of increasing diameter, with flared transitions rather than 
vertical plates wherever possible between two sections of pipe (the pipe bent 
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out at a 30° angle to meet the adjoining segment), and by a thin window at 
z = 120 in through which all collision primaries entering the detector would 
pass. The steps in the Al pipe are given in Table \i.2[ and the beampipe and 
abort are pictured in Fig. p.2| . The new pipe was installed in February 1995, 
and was in place for the majority of the data collected. 

The main ring, which was used for accelerating protons which were then 
used in production of the anti-protons, and a main ring abort pipe were 
present in the CO collision region, both approximately 20 in above the Teva- 
tron pipe. These pipes were steel, on the order of 1/4 in thick. The main ring 
and abort were about 6 in and 4 in in diameter, respectively. The accelerat- 
ing protons in the main ring interacted with particles present in gas inside 
the pipe, producing blinding flashes of background in the detector, so that 
the triggers had to be gated off when the main ring protons passed through 
the collision region. In events which were triggered on, such beam-gas in- 
teractions also created hits in the upstream scintillator used for diffractive 
tags (which will be discussed in Sec. |2.2|) . Interactions of collision primaries 
with the main ring pipes produced background tracks in the detector. And 
accelerator studies using the main ring produced large bursts of radiation 
which badly damaged some of the electronics. This was the greatest source 
of radiation damage because the spread in momentum of beam particles in 
the (older) main ring was much greater than that for the Tevatron beams, 
and those particles with less-than-beam momentum often interacted in the 
pipe or escaped the pipe and interacted elsewhere. 



2.2 The MiniMax detector 



2.2.1 MWPC telescope 

Particles were tracked through a telescope of 24 multi-wire proportional 
chambers on the downstream proton side of the collision point. The ac- 
ceptance of these chambers was approximately a circle of radius 0.6 units 
in 7] — (j) space centered on rj ^ 4 (see Fig. |2.3|). The chambers were de- 



signed for past experiments at CWRU [34]. Figures 2.4 and 2.5 show dia- 



grams of the wire chamber design. Each chamber had an active volume of 
12.8 in X 12.8 in x 0.375 in containing 128 parallel wires with a spacing of 
0.1 in. The anode wires were gold-plated tungsten with a 0.0008 in diame- 
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ter, and the cathode plane consisted of an evaporated Al film on one side of 
0.001 in-thick mylar. The gap was 3/16 in. Copper clad surfaces 1/16 in from 
the signal plane served as guard rings, and 3/16 in from the signal plane as 
the contact for the high-voltage plane. The 80% Ar, 20% CO2 gas used was 
circulated freely on both sides of the aluminized mylar. The outer faces were 
sealed with 0.005 in-thick mylar. After tests in early 1995 determined that 
some of the chambers were sensitive to light, 0.002 in-thick black plastic film 
(kevlar, or similar material) was added to the front and back faces of the 
chambers. 

The chambers were held together by l/8in- and l/16in-thick frames made 
of G-10 which were epoxied together. Grooves 0.004 in-wide were etched in 
the printed circuit board face of the signal plane to keep the signal wires 
aligned during the epoxy process. The 1/16 in-thick G-10 circuit board ex- 
tended 1.85in below the frame, which was lin in width, so that the completed 
chambers had a width of 14.8 in, and a height of 16.65 in. 

Two types of readout electronics were used. Half of the chambers used 
nano metrics supplied by Fermilab, which did not give pulse- height informa- 
tion, but were very reliable. The other half were equipped with cards made 
by the University of Michigan. These had 12 bits to digitize the signal from 
each wire, and gave pulse-height information about the amount of charge 
deposited on the wire. A certain chip in the amplifier cards was had to be 
replaced frequently due to radiation damage, which often led to runs which 
did not have all MWPC's functioning .0 

The chambers were held in position by aluminum stands which allowed 
them to be rotated at various angles. The original stands each held four 
chambers. New stands were made in September of 1995 to hold eight cham- 
bers in order to compress them closer together. An end-on view of a chamber 
stand is pictured in Fig. p?6| , along with other pieces of apparatus. 

Chamber orientations for various periods of running are given in Tables 
2^2. 6|. The alignments during the first year or so of running utilized large- 



angle stereo for three-dimensional resolution. Each chamber was rotated by 
a different angle in order to reduce the number of potential reconstructed 



tracks as described below. This philosophy was reconsidered [35], and for 



^For a brief period in Dec. 1995, the Tevatron ran at y/s = 630 GcV. However, all 
of the data collected by MiniMax at that energy did not have information from the wire 
chambers which were read out using the Michigan electronics because these were under 
repair at the time. 
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reasons discussed below, the chambers were reconfigured in February 1995. 

The final orientation of the chambers was defined in a coordinate system 
where, looking in the positive z direction, the (m, f )-axes were a 45° counter- 
clockwise rotation of the (x, ?/)-axes about the z-axis. [We also redefine the 
azimuthal angle as = tan"^ (f/w).] In the chamber configuration used for 
the majority of the data taking, three of the front chambers and eight of 
the rear had their wires aligned perpendicular to the w-axis ("m chambers"), 
three of the front chambers had wires aligned within 15° of normal to the v- 
axis ("f chambers"), and the remainder were rotated by small angles (4-15°) 
from the u-chamber orientation ("u' chambers"). 

A configuration with wires in half the chambers aligned in a given direc- 
tion and the other half perpendicular gives equal resolution in both directions. 
However, the number of potential tracks considered by a track-finding algo- 
rithm goes like the square of the number of real tracks. For example, two real 
tracks will produce two hits in a front chamber and two in a back chamber, 
leading to four potential track candidates if such candidates are found using 
only the two chambers. In a situation like that of MiniMax, where there 
are many random hits due to pipe shower, confirming hits in other chambers 
can make these extra track candidates look like actual tracks. The opposite 
extreme is a chamber configuration with all wires parallel. This does not 
give any information about the location of a track in the direction along the 
wires, but greatly reduces the number of potential fake tracks. The other 
great advantage to parallel wires, as opposed to large-angle stereo at random 
angles, is that tracks can be reconstructed by eye in a two-dimensional dis- 
play of the hit wires, which can be quite useful for studying the properties 
of, e.g., tracks coming from interactions in the beampipe. 

The MiniMax Collaboration determined that the configuration discussed 
above, employing small-angle stereo, was a good compromise between the 
two extremes because it allowed for some resolution in both directions, the 
ability to view wire hits in about half the chambers, and reduced probability 
of finding fake tracks. Note also that the front chambers had three v chambers 
for increased resolution of charged track^ and of their intersection with the 
z-axis (i.e. the collision point); the number of wire hits, especially of those 

charged track is defined as the reconstructed track from a charged particle which 
appears to originate in the collision. 
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due to pipe shower, was much smaller in the front chambers than in the back, 
so that the reduction of fakes by small-angle stereo would not have been as 
much of an advantage as it was in the rear chambers. 

Another major reconfiguration was performed in the fall of 1995 in which 
the chambers were moved closer to the collision point. The reason was that 
pipe shower and other sources of background increased towards the rear of 



the detector. This will be discussed further in Sec. 3.2.1 



2.2.2 Converter 

In order to detect photons, converter was inserted behind the eighth chamber. 
During early running (through 1994), the converter was a stationary plane of 
lead, 6inx4inx lXo,0 centered on (5.25in, 6.20in, 198in), level with the floor. 
For runs after Feb. 1995, the movable converter was an 8 in x 8 in square, 
rotated by 45° relative to the floor, so that its edges were parallel to the {u, v) 
axes. Various thicknesses and materials were used to study systematics in 
detecting photon conversions; these were: 1, 2, 0.5 Xq lead, 1 Xq iron, and 
no converter (converter in the "out" position). For runs before chamber 
compression, the lead was located at (7.7 in, 0.1 in, 184.8 in) for "lead-in" 
running, and approximately 20 in farther from the beampipe in x for "lead- 
out" runs. For the production runs, 1 Xq lead was used about equally often 
in the lead-in position (5.13 in, 5.13 in, 150 in), and in the lead-out position. 



2.2.3 Trigger scintillator 

Scintillator counters were used to trigger on particles passing through the 
detector, signaling a collision. The counters had both ADC (analog to dig- 
ital converter) readout, which gave information about the energy deposited 
by particles passing through the scintillator, and TDC (time to digital con- 
verter) information about the time-of flight of the particles. The signature of 
beam-beam (rather than beam-gas) collisions is correctly-timed hits in the 
scintillator arrays. 

■^One radiation length (1 Xq) is defined as the mean distance over which a high-energy 
electron loses all but 1/e of its energy by bremsstrahlung, and varies with atomic number 
and mass of a medium. Also, for very high-energy photons, the e'''e~ pair-production cross 
section is given by cr « ^(A/XqNa), where A is the atomic mass and Na is Avagadro's 
number p^ . 
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The scintillator counters which were interspersed among the chambers 
were referred to as A-E, or collectively as "alphabet counters" . The B and D 
counters which played a role in the usual trigger were located one directly be- 
hind the lead and one directly behind the last chamber. Each was comprised 
of two 8 in X 16 in X 0.5 in pieces which together formed a square, which was 
rotated into the (u, v) frame in a similar manner as the lead. The C counter 
was of the same type as B and D. In early runs, B and C were located on ei- 
ther side of the converter, and signals from both were required for the trigger, 
significantly reducing the probability of triggering on an event with a photon 
conversion but no charged tracks. The smaller A counter, located in front of 
the MWPC telescope, and the larger E counter behind the calorimeter were 
also sometimes used in the trigger. The history of configurations is given in 
Tables 

The trigger also required a hit in scintillator counters on the downstream 
anti-proton side of the collision ("pbar counters"). This array contained four 
12.5 in X 2 in X 1 in pieces and four 12.5 in x 4 in x 1 in pieces arranged to 
form an 18.5 in square with a 6.5 in square hole in the middle which was 
centered around the 6 in-diameter beampipe at 2; = —81 in, as shown in Fig 
|2.6| . Early running required a hit in an identical array on the opposite side 
of the collision at z = 83 in ("p counters"), but the information provided by 
these and by the alphabet counters was redundant, and the p counters were 
removed when the new beampipe was installed. 

Histograms of the ADC values for some of the counters are shown in Fig. 
\2.T\ . The trigger was designed to record events in which minimum-ionizing 
particles passed through the B, D, and pbar counters. [Any moderately- 
relativistic charged particle other than an electron is a minimum-ionizing 
particle (mip), and loses energy in a medium mainly through ionization. 
The amount of energy lost depends on the thickness of the material and 
the velocity of the particle; the mean energy lost over a given distance is 
described by the Bethe-Bloch equation [BBI.] Plots of the pbar ADC's show 



these mip peaks, which indicate the number of minimum-ionizing particles 
which passed through. These peaks are not seen clearly in the alphabet 
ADC's because of degradation of the scintillator, causing poorer correlation 
between the amount of light created by the mip track and that collected by 
the photomultiplier tubes; however, by increasing the voltage on the counters, 
triggering on minimum-ionizing particles was still possible. An entry in the 
log book from January 7, 1996 notes this, along with particular mention of 
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the poor performance of the top D counter. 

A signal from the accelerator division reported the occurrence of a beam 
crossing at CO. At this signal, counting began in the TDC's, and the time 
when a charged particle was detected in the scintillator was recorded. If no 
charged particle was detected after 2^^ counts, a hit was recorded in the last 
bin. Plots of the TDC values for the B, D, and Pbar counters (summed over 
all pieces of scintillator) are shown in Fig. 2^ . One TDC count is equivalent 
to 0.25 ns for all counters. However, the offsets for each counter are different 
due to differing lengths of cable, so that a particular TDC value does not 
correspond to the same amount of real time for all counters. Also shown in 
Fig. |2]^ is the TDC for the bottom D counter plotted against the TDC for 
one of the pbar counters. The dark region indicates the timing for beam- 
beam collisions. Hits which occur earlier or later in one of the counters are 
due to something other than collision primaries. For example, a beam-gas 
interaction involving a proton and occurring before the proton bunch reaches 
the collision region could send particles into the pbar counters at times earlier 
than those for beam-beam collisions. Also, secondary interactions in the 
detector could cause hits at later times. 

Since the counters on the detector side were behind the lead converter, 
we were able to trigger on events with a photon conversion regardless of 
the presence of charged tracks. In fact, the large frequency of interactions 
with the beampipe which sent particles into the detector resulted in many 
triggered collisions which produced no primary charged particles or photons 
in the acceptance. 

The trigger rate during normal running conditions (DO luminosity of 10^° 
cm~^s~^) was around 300 Hz. Only a fraction of the events which were 
triggered on were recorded; events were written at a rate of approximately 
30 Hz. The fraction of triggers not due to beam-beam collisions was less than 
about 5%, and was mainly due to proton-induced beam-gas interactions. 

Mean trigger rates and Tevatron running conditions for some of the Jan- 
uary 1996 runs are given in Table |2.10| . The BDpbar delay is an indirect 
measure of the time difference of signals from pbar and alphabet counters, 
and different values lead to triggers on beam-beam or beam-gas collisions. 
The p and pbar currents are the numbers of protons and anti-protons cir- 
culating in the ring. The raw trigger rate requires ADC counts from the 
B, D, and pbar counters which are high enough to indicate the presence of 
collision-produced particles, and the beam trigger rate additionally requires 
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coincidental signals from the TDC's of these counters. (The alphabet rate is 
not given because it was not recorded correctly in these runs.) The ratio of 
beam trigger rate to DO luminosity seems to vary more in the low luminosity 
runs than in the runs with higher luminosity. The ratio of beam to raw rates 
is a measure of the amount of raw triggers due to beam-gas interactions. For 
the last few runs in Table |2.1CI| , the beam rate is reported as being higher 
than the raw rate, which is clearly impossible. The log book notes a "dou- 
ble pulsing" in the trigger rate beginning in run 1128, and during run 1137, 
"beam trig / raw trig is ramping to 2.0 then back down over ~ 5 min period 
- has occurred 4 times over last 200k events" . The source of the problem is 
unclear. 

The cross section seen by the MiniMax detector was estimated using the 
cross sections determined by CDF. These are ||3^ 80.03 ± 2.24 mb total cross 
section, 19.70 ± 0.85 mb elastic, and 9.46 ± 0.44 mb single diffractive (this is 
the sum of cross sections for diffractive dissociation of the proton and of the 
anti-proton). The non-single diffractive inelastic cross section is therefore 
50.87 mb. Some single diffractive events (most likely with high diffractive 
mass) may have been triggered on, and not all double diffractive events pass 
the trigger. The cross section of events which actually passed the trigger was 
determined from the trigger rate and luminosity. The luminosity at CO was 
taken as that at DO corrected for differences in the magnetic architecture in 
the two regions and the fact that the bunches which collide at CO are not the 
same pairs as those which collide at DO. The parameter /? is a measure of 
how tightly focussed the beams are at a given point. The ratio of the /3's at 
CO and at DO is approximately 0.35/72, and is further corrected by a factor 
of 1/0.7 to take into account the variation of (3 over the luminous region at 
DO. The ratio of trigger rate to DO luminosity was typically around 0.30 mb. 
Therefore we estimate that the observed cross section was approximately 
43 mb. 



2.2.4 Electromagnetic calorimeter 

A 28-module electromagnetic calorimeter was located behind the MWPC's, 
with the face aX z = 282in, centered on x = 8.49in, y = 7.71 in. Twelve of the 
cells were 3.8 x 3.8 in^ with 30 layers each of 4.88 mm-thick lead and 6 mm- 
thick scintillator, which combined for a total length of about 27 Xq. These 
were made by the University of Michigan |^ and were used in Fermilab 
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El-756. (Before 1995, the calorimeter consisted of 16 Michigan cells with 
the same total cross sectional area, and with the face at z = 288 in.) Four 
10cm X 10cm X 30cm cells each had 4 photomultiplier tubes, which effectively 
made 16 5 cm x 5 cm x 30 cm modules. The cells, which were cut from longer 
hadronic calorimeter cells made by Wayne State University, were lead with 
47x47 scintillating fibers of diameter 1mm spaced 0.213cm apart and running 
the length of the module The effective radiation length was 0.78 cm, 

so that the modules were about 38.5 Xq long. The calorimeter was placed 
on a stand so that it also was rotated by 44.5°. The smaller cells occupied 
positions close to the beampipe in order to cover roughly the same area in 
T) — (p space as the larger cells, as seen in Fig. |2.9| (recall that larger values of 
T) are closer to the beampipe at a given z) . This is desirable since the density 
of particles is roughly uniform in lego space. 



2.2.5 Upstream tags 

In order to tag diffractive events, scintillator and hadron calorimetry was 
placed upstream of the collision point. A view of the upstream region is shown 
in Fig. |2.10| . At z ~ —25 m, two hadron calorimeter modules, similar to the 



smaller cells of the electromagnetic calorimeter except having dimensions of 
lOcmx lOcmx 117cm, were positioned to detect leading anti-protons and anti- 
neutrons. The machine magnets bent the p's with less-than-beam momentum 
by greater angles than the beam, and the n's not at all. The anti-protons 
detected here had a Feynman x of Xi? ~ 0.5, i.e. their momentum is about 
half that of the beam particles. Farther upstream at about z = —60 m, four 
scintillator counters were placed to detect showers from anti-protons with 
xp ^ 0.9 that interact in what are called the kicker magnets. If at least one 
of these scintillator is hit, the event has a "ktag". Note that these counters 
did not only see diffractive anti-protons, but also products from beam-gas 
interactions in the upstream region. 



2.3 Run history 



The concept of a detector with good coverage in the forward region which 
ultimately led to the MiniMax experiment originated in an initiative for 
a full-acceptance detector (FAD) at the SSC. A proposal was later 
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submitted for a maximum-acceptance experiment at Fermilab (MAX) |^ 
which was designed to investigate some of the physics goals of FAD. A small 
test /experiment with a detector covering only the forward region (MiniMax, 
Fermilab T-864) was approved in the spring of 1993. From late 1993 to early 
1995, the MiniMax detector was installed, and collider data was taken with 
8, 12, 16, and finally 24 chambers in large-angle stereo. My work on the 
experiment began in September 1994. 

The new beampipe was installed in February of 1995, and when the cham- 
bers were re-installed, they were given the new small-angle-stereo orientation 
with wires in 11 of the 24 chambers aligned parallel to each other for rea- 
sons discussed in the previous section. During the following five months, 
short opportunistic runs were taken with various thicknesses of lead and iron 
converter, and two different window thicknesses, in which a total of about 
3.5 X 10^ events were collected. The diffractive tags were commissioned in 
May and were available in more than 10^ events. The chambers were again 
removed and reconfigured in the fall of 1995 in order to compress them closer 
to the interaction point, away from some of the blinding pipe showers farther 
downstream, as will be discussed in Sec. |3.2.1| . 

Production running in January of 1996 yielded 4.2 x 10^ events, 2.7 x 10^ 
of which had IXq lead in, and the remaining 1.5 x 10^ events were run with the 
lead out. The luminosity in these runs was lower than in most earlier runs by 
about 1-2 orders of magnitude. This also led to fewer beam-gas interactions 
which, among other effects, meant that the sample of events with a ktag had 
a much higher fraction of diffractive events. Some information about the 



runs used in the analysis described in this work is given in Table p. 10 



The MiniMax detector was decommissioned the following spring. 
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Figure 2.1: The MiniMax detector before compression of the MWPC tele- 
scope. 
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I. 



Figure 2.2: New Tevatron beampipe and abort at CO and enlarged view of 
the region from z ~ 110 — 190 in. 
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Figure 2.3: Acceptance in lego space. 
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Figure 2.4: MWPC front view. 
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Figure 2.5: MWPC top view. 
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Figure 2.6: Pieces of apparatus including pbar and alphabet scintillator, 
chamber stands, and the original calorimeter. 



27 




Figure 2.7: ADC spectra from some of the trigger counters. 
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Figure 2.8: TDC spectra from some of the trigger counters. 




Figure 2.9: Electromagnetic calorimeter cells in lego space; the vertical axis 
is (f) and the horizontal axis is 77. 
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Figure 2.10: A view of the detector including the upstream region. 
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Table 2.1: Old Tevatron beampipe and abort at CO. 
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Table 2.2: New Tevatron beampipe and abort at CO. 
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6 


1.96328 


10.8733 


-0.442493 


208.952 


Michigan 


7 


-3.03528 


14.9853 


6.86178 


212.912 


Michigan 


8 


1.41509 


7.64913 


1.50400 


216.934 


Michigan 


9 


-0.400473 


4.92871 


10.5747 


247.129 


Michigan 


10 


-2.21307 


14.7361 


13.0676 


250.841 


Michigan 


11 


2.23920 


14.9585 


3.01641 


256.658 


Michigan 


12 


-2.74165 


16.9554 


10.4961 


260.679 


Michigan 



Table 2.3: Chamber alignment from 2/16/94. The location and orientation 
of each chamber is defined by the {xq, i/o, Zq) coordinates of the midpoint of 
wire number zero, together with the angle with respect to the x-axis of the 
vector pointing in the direction of increasing wire number. 
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nanometrics 
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Michigan 
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207.02 


Michigan 
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2.0755 


10.282 


1.746 
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Michigan 
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Michigan 
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nanometrics 
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nanometrics 
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0.2893 
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nanometrics 


120 


-1.2743 


6.283 


14.190 


241.58 


nanometrics 


121 


-1.4826 


7.788 


14.272 


246.58 


Michigan 


122 


-3.0389 


14.772 


8.508 


251.89 


Michigan 


123 


-0.7090 


3.398 


12.185 


254.58 


Michigan 


124 


-2.2739 


12.464 


12.906 


256.14 


Michigan 



Table 2.4: Chamber alignment from 11/23/94. 
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Michigan 


113 


-2.3558 


11.167 


11.409 


212 


84 


nanometrics 


114 


-2.5271 


12.242 


10.461 


216 


84 


nanometrics 


115 


-2.3580 


12.240 


11.516 


221 


09 


nanometrics 


116 


0.9668 


3.447 


1.569 


226 


66 


nanometrics 


117 


-2.3592 


11.679 


11.749 


237 


72 


nanometrics 


118 


-2.2361 


11.216 


12.239 


241 


72 


nanometrics 


119 


0.7847 


2.616 


2.686 


245 


84 


nanometrics 


120 


-2.4771 


12.387 


11.188 


249 


72 


nanometrics 


121 


-2.2556 


11.497 


12.241 


256 


15 


Michigan 


122 


0.7104 


2.377 


3.154 


260 


15 


Michigan 


123 


0.7801 


2.779 


2.820 


264 


40 


Michigan 


124 


0.7818 


2.775 


2.816 


265 


97 


Michigan 



Table 2.5: Chamber alignment with 11 chambers parallel from 3/19/95. 
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-2.3370 


10.970 


11.702 


172.84 


Michigan 


112 


0.5794 


1.535 


3.508 


176.09 


Michigan 


113 


-2.3370 


11.199 


11.781 


178.84 


nanometrics 


114 


-2.5063 


12.360 


10.771 


181.84 


nanometrics 


115 


-2.3405 


11.424 


11.812 


185.09 


nanometrics 


116 


0.9669 


3.366 


1.695 


189.65 


nanometrics 


117 


-2.3370 


11.602 


12.036 


194.29 


nanometrics 


118 


-2.2201 


11.327 


12.524 


197.29 


nanometrics 


119 


0.8029 


2.814 


2.963 


200.29 


nanometrics 


120 


-2.4574 


12.648 


11.415 


203.29 


nanometrics 


121 


-2.3527 


12.056 


12.114 


206.16 


Michigan 


122 


1.0420 


3.560 


1.613 


209.16 


Michigan 


123 


0.7837 


3.140 


3.048 


212.41 


Michigan 


124 


0.6266 


2.681 


3.748 


215.41 


Michigan 



Table 2.6: Chamber alignment after compression from 
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counter 


X 


y 


z 


height 


width 


angle 


A 


-2.8 


-2.2 


108. 


5.5 


5.5 


0.26 


Bl 


-5.9 


-3.7 


155. 


8.0 


16.0 


0.26 


B2 


-3.8 


3.9 


155. 


8.0 


16.0 


0.26 


CI 


-5.9 


-3.7 


223. 


8.0 


16.0 


0.26 


C2 


-3.8 


3.9 


223. 


8.0 


16.0 


0.26 


Dl 


-5.9 


-3.7 


269. 


8.0 


16.0 


0.26 


D2 


-3.8 


3.9 


269. 


8.0 


16.0 


0.26 


E 


-12.2 


-4.0 


322. 


13.0 


13.0 


1.04 



Table 2.7: Configuration of alphabet counters for runs with the old beampipe, 
through 1994. 



counter 


X 


y 


z 


angle 


A 










Bl 


5.1 


10.8 


189. 


45. 


B2 


10.8 


5.1 


189. 


45. 


CI 


5.0 


10.7 


180. 


45. 


C2 


10.7 


5.0 


180. 


45. 


Dl 


6.7 


11.3 


280. 


41.7 


D2 


12.0 


5.3 


280. 


41.7 


E 


-12.2 


-4.0 


322. 


0. 



Table 2.8: Configuration of alphabet counters during running with the new 
beampipe from 2/95-7/95. 



counter 


X 


y 


z 


angle 


A 


2.8 


2.8 


264. 


45. 


Bl 


6.1 


11.8 


279. 


45. 


B2 


11.8 


6.1 


279. 


45. 


CI 


6.0 


11.7 


270. 


45. 


C2 


11.7 


6.0 


270. 


45. 


Dl 


7.2 


7.6 


157. 


27.4 


D2 


10.9 


0.5 


157. 


27.4 


E 


-12.2 


-1.0 


322. 


0.0 



Table 2.9: Configuration of alphabet counters during runs with compressed 
MWPC telescope. 
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run 


date 


lead 


number 


BDpbar 


P 


pbar 


DO 


raw 


beam 


beam/ 


beam/ 


number 






events 


delay 


current 


current 


luminosity 


rate 


rate 


DOlum 


raw 










(ns) 


(10^) 


(109) 


(103°cm-2s-i) 


(Hz) 


(Hz) 


(mb) 




1089 


1/18 


out 


50112 


8 


105. 


83. 


0.91 


269. 


258. 


0.28 


0.96 


1093 


1/18 


out 


103061 


8 


95. 


79. 


0.63 


192. 


190. 


0.30 


0.99 


1096 


1/18 


in 


249986 


8 


38. 


30. 


0.14 


43. 


42. 


0.30 


0.97 


1099 


1/18 


in 


174775 


8 


37. 


29. 


0.09 


32. 


32. 


0.33 


0.99 


1103 


1/19 


in 


67117 


8 


34. 


27. 


0.05 


18. 


18. 


0.35 


0.98 


1104 


1/19 


out 


50461 


8 


34. 


26. 


0.04 


13. 


13. 


0.36 


0.98 


1108 


1/19 


out 


249990 


8 


37. 


30. 


0.15 


45. 


45. 


0.31 


0.99 


1109 


1/19 


in 


250081 


8 


36. 


29. 


0.11 


34. 


34. 


0.32 


0.99 


1110 


1/20 


in 


249967 


8 


35. 


28. 


0.08 


27. 


26. 


0.32 


0.98 


nil 


1/20 


out 


234361 


8 


34. 


27. 


0.06 


20. 


19. 


0.33 


0.98 


1123 


1/21 


out 


37126 


8 


61. 


74. 


0.45 


128. 


128. 


0.29 


1.00 


1124 


1/21 


out 


91958 


8 


21. 


14. 


0.05 


13. 


13. 


0.28 


1.00 


1125 


1/21 


in 


251434 


8 


20. 


13. 


0.03 


9. 


9. 


0.36 


1.00 


1126 


1/22 


in 


47944 


8 


19. 


13. 


0.02 


8. 


7. 


0.35 


1.00 


1127 


1/22 


in 


28933 


6 


19. 


13. 


0.01 


3. 


3. 


0.38 


0.99 


1129 


1/22 


out 


53767 


10 


70. 


85. 


0.55 


195. 


204. 


0.37 


1.05 


1132 


1/22 


out 


52402 


4 


68. 


80. 


0.41 


138. 


146. 


0.36 


1.06 


1137 


1/22 


in 


249977 


8 


65. 


74. 


0.26 


128. 


137. 


0.52 


1.07 


1139 


1/23 


in 


153014 


8 


60. 


67. 


0.16 


75. 


74. 


0.46 


1.00 



Table 2.10: Running conditions from some January 1996 runs. 



Chapter 3 
Simulations 



3.1 PYTHIA 

Minimum-bias collisions at ^/s = 1.8 TeV are simulated using PYTHIA ver- 



sion 5.702 and JETSET 7.401 |42, 43]. The combination of PYTHIA and 



JETSET provide a commonly-used event generator for high-energy collisions 
of elementary particles based on the parton model |Q, and will be hereafter 



referred to as "PYTHIA" . The physics of the collisions that PYTHIA deals 
with includes 

the parton distribution functions of the beam particles, i.e., the flavor of 
the quarks and the fraction of energy carried by the constituent quarks 
and gluons of the incident proton and anti-proton, 

any radiation from the beam particles, such as the emission of a gluon 
from a quark, which leads to an initial-state shower, 

the hard process between the incoming partons from each shower which 
produces outgoing particles, 

any final-state showers produced by outgoing partons, 

the beam remnants, which must return to color-neutral states after 
losing the interacting partons, 

the hadronization of outgoing partons - confinement requires that they 
form color-neutral hadrons, 



39 



40 



and the decay of any unstable outgoing particles. 

For the hadronization, PYTHIA uses a string fragmentation model, specif- 
ically what is called the Lund model PSI. Its basis is a string connecting a 



quark and an anti-quark which stretches as the partons move away from each 
other until it breaks by creating a new quark-anti-quark pair. This contin- 
ues until only on-mass-shell hadrons remain, where a hadron is a color-singlet 
pair connected by a small piece of string. 

Default values are taken for all parameters except that particles with a 
mean invariant lifetime cr greater than 1 cm are not decayed. This allows 
the decays of K° and A° particles to be studied later. 

For non-single diffractive inelastic pp collisions at 1.8 TeV, PYTHIA gives 
mean numbers of particles into the acceptance per event of 0.61 vr^'s, 0.15 
other charged particles, and 0.70 7's. The inclusive pseudorapidity distri- 
butions, dN/di], for charged pions, all charged particles, and photons are 
shown in Fig. |3T^. (An interesting note is that the larger fraction of particles 
with 6 < ?7 < 10 in the charged distribution relative to the charged-pion dis- 
tribution is due to leading protons and anti-protons from double diffractive 
processes.) 

The cross sections for various types of events included in the minimum- 
bias events generated by PYTHIA are given in Table Note that compared 
to the cross sections reported by CDF (Sec. p.2| ), PYTHIA underestimates 
the total and elastic cross sections, and overestimates the single diffractive. 



3.2 GEANT 

The particles generated in a collision are then taken as input into a GEANT 
simulation (GEANT version 3.21 [^) where they are propagated through the 
material of the detector. The initial position of the collision primary particles 
is taken to be m = f = 0, and z given by a Gaussian distribution with a mean 
of 7 in and standard deviation 16.2 in, in order to reproduce the variance 
of the collision point seen in the data. Cross sections and simulations of 
processes such as hadronic interactions, electromagnetic processes, ionization 
by charged particles, multiple scattering, and decays in flight are included 
in GEANT. Default routines are used for all such processes, as well as the 
default energies below which particles are not tracked (IMeV for photons and 
electrons, lOMeV for hadrons and muons). The energy deposited in material 
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due to interactions is given at each step, and the secondary particles produced 
are added to the hst of particles to be tracked through the detector. 

The energy deposited in the detector elements is used to model the re- 
sponse to particle interactions seen in the real data. At least 4.4(2.2) MeV is 
required to be deposited in the 1 in (0.5 in)-thick scintillator trigger counters 
in order to keep the event. Figure |3]^a shows the energy deposited in the 
1 in-thick scintillator counters for minimum bias events, and Fig. |3.3| a shows 
the energy deposited by single charged pions (minimum-ionizing tracks). In 
Fig. |3.2| a, the peak for one minimum- ionizing particle passing through the 
scintillator can be seen around 4.5 MeV, and the smaller peak for two mip's 
at 9 — 10 MeV. The first minimum-ionizing peak starts at about 4 MeV, and 
has a maximum at 4.5 MeV. The low energy background is fairly insignifi- 
cant. Therefore, any cutoff energy between these values is reasonable. The 
value 4.4 MeV is chosen to be as large as possible in order to increase the 
mean number of wires hit in the MWPC's ("NHITS") for triggered events 
without missing a significant part of the mip peak. (The reason for this will 
be given in Sec. |3.2.1| .) 

For the GEANT simulation, the gas in each MWPC is segmented into 
128 pieces corresponding to the volume which surrounds each wire in the 
actual chambers. The wires are not included in the simulation. Plots similar 
to those for the scintillator energy are given for the chamber pulse heights 
(energy deposited in a segment) in Fig. ^]2|b and Fig ^]3|b. For single charged 
pions, the low energy background which falls off around 0.3 — 0.4 keV where 
the mip peak begins is presumably due to tracks which pass through more 
than one segment and deposit some fraction of the minimum- ionizing energy 
in each segment. For minimum-bias events, the low energy background is 
much higher and extends into the minimum-ionizing peak. The cutoff value 
chosen to represent a hit wire is 0.4 keV, where the background falls to a 
level comparable to the mip contribution. Note that in the actual MWPC's, 
ionized electrons created by the charged tracks drift to the positively charged 
wires and the amount of charge collected on a wire determines the pulse 
height. This may not be linearly related to the energy deposited in the gas; 
however, direct comparisons of pulse heights in the GEANT and in the data 
discussed in Sec. |3.2.1| show that this is not an unreasonable assumption. 

The energy deposited in each cell of the lead-scintillator electromagnetic 
calorimeter is recorded, so that it can be used to calibrate the signals from 
the actual calorimeter. 
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The other crucial detector element which is included is converter of vari- 
ous thicknesses. 

Because of the large observed background of particles due to interactions 
in material in the forward (detector) region, many objects not related to the 
detector are included in the simulation. These include the Tevatron beam 
pipe and abort pipe, the main ring and its abort, support stands, vacuum 
pumps, the concrete floor, etc. In spite of all these additions, the mean 
number of wires hit in each MWPC is lower than that in the real data by 
about a factor of two (see Figs. ^]4|, ^l5|) . 



3.2.1 GEANT NHITS study 

The NHITS distribution is dependent on the pulse-height and trigger cuts. 
A lower pulse-height cut allows lower-energy hits to be counted as hits by 
actual tracks, thus increasing the NHITS in each event. A higher cut on 
energy deposited in the scintillator counters in order to trigger (i.e. to keep 
that event) tends to cut out events with less energy or number of tracks going 
into the detector, so that the mean NHITS of the remaining events is higher 
than that for all events. The cuts made on the GEANT data are intended to 
match those in the real data, namely to make cuts based on the signatures 
of real charged tracks going through the chambers and scintillator. However, 
we tried varying these cuts through all reasonable (and even unreasonable) 
values to see if the NHITS distribution of the data could be reproduced by 
the GEANT simulation. Only 12 of the 24 MWPC's read out pulse-height 
information, so we use only those chambers in the following analysis. The 
pulse height is calibrated in terms of charge deposited on the wires, not 
energy deposited in the gas as it is in GEANT, and is not the same in all 
chambers. The first step, therefore, is to adjust the energy scales of the 
chambers in the GEANT so that the mean of a Gaussian fit to the peaks 
of both the real and GEANT pulse- height distributions are equivalent. The 
result is shown in Figs. |3.6| , |3.7| . The mean NHITS in these 12 chambers 
is 93 for run 867. It is clear in Table 13.21 that there is no combination of 



pulse-height and scintillator cuts which gives a mean NHITS this high in the 
GEANT data. For the lowest possible pulse- height cut (pulse height > 0), 
and the highest shown cut on energy deposited in the scintillator (8 MeV, 
which practically requires two charged tracks rather than one), the mean 
NHITS in those 12 chambers is only 79. 
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The next attempts to find an explanation for tlie lower NHITS in GEANT 
include the following: 

1. Increasing the thickness of the beampipe from 0.03 in to 0.035 in to 
check for GEANT problems with the boundary; 

2. lowering the energy thresholds of particles for GEANT to track to 1/2 
the default value; 

3. changing the defaults of GEANT so that delta rays are produced; 

4. supposing that the collision point is not where we think it is relative 
to the pipe, and therefore changing the z of the collision point to 10 in 
closer to the pipe, or changing x and y by 1 in each towards the pipe. 

None of these have much effect on the NHITS. Next we tried adding the main 
ring and correcting the description of the abort pipe to include the change of 
material to steel at z > 275 in, which increased the mean NHITS by about 
20%. Encouraged by this, we added more of the material nearby the detector: 
2ft-thick concrete floor, tunnel wall, chamber stands, Tevatron-beampipe and 
main-ring support stands, and vacuum ion pumps. 

The effect on NHITS of material close to the wire chambers was studied 
in more detail by removing that material. The beampipe is the largest source 
of non-collision-point tracks hitting the chambers (Fig. The abort pipe 

and the G-10 frames of the wire chambers were also studied (Figs. p.9| , |3.10|) . 



The rear chambers are hit by tracks from these sources more than the front; 
in fact, when all three of the sources are removed, the mean number of hits 
per event in each chamber is approximately a constant between 3.5 and 4 
(Fig. p. lip , as opposed to the roughly linear increase in hits going from 



about 4 in the front chamber to more than 8 in the last, as in Fig. 
(This is, of course, for lead-out runs. For lead-in runs, the number of hits 
in the chambers behind the lead is much greater than that in the front due 
to conversions of primary photons.) In the real data, the number of hits 
per chamber increases faster than linearly towards the rear of the detector. 
The addition of the beampipe supports in the GEANT simulation reproduces 
this effect fairly well. These observations are what led to the compression of 
the wire chambers which moved them closer to the collision point before the 
production running. 
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It has been noted that the multiphcity distributions of \ow-pT charged 
particles in PYTHIA were incorrectly extrapolated from collider data at high 
Pt- Since the MiniMax detector is more sensitive than conventional collider 
detectors to low-p^ particles, these can potentially contribute significantly to 
the NHITS. We tried using HERWIG, another commonly-used event genera- 
tor (especially for jet studies) to see if the NHITS distribution might be better 
reproduced. Comparisons of dN/ drj from non-diffractive PYTHIA events and 
HERWIG events show that the mean number of particles produced is higher 
for PYTHIA (see Fig. |]T|for dN^h/drj). The NHITS distribution obtained 
using HERWIG was not significantly different. 

Another attempt to increase the mean NHITS in the simulations was to 
change the defaults for multiple interactions in PYTHIA. The occurrence of 
hard interactions of more than one parton pair in a hadronic collision is not 
well understood, and PYTHIA provides several models. The default is that 
the probability of multiple interactions is equal for all events, with a sharp 
Vi^min cut-off. It was suggcstcd to us that the option with "multiple interac- 
tions assuming a varying impact parameter and a hadronic matter overlap 
consistent with a double Gaussian matter distribution . . . with a continuous 
turn-off of the cross section at p^o" (Ref. [0, p. 222) would give larger mul- 
tiplicity fluctuations, and therefore possibly larger mean NHITS. The effect 
of using this model is also negligible. 

A study was done by the ALICE Collaboration on MWPC's which are 



intended to be used in a muon spectrometer at the LHC [^. The cham- 



bers are remarkably similar to the ones used by MiniMax, and therefore the 
results obtained in the study may be relevant to the performance of the Min- 
iMax chambers. A pion beam was aimed at a lead absorber, and background 
particles from interactions in the lead were detected in a wire chamber lo- 
cated at the position which was predicted by simulations as the location of 
the shower maximum. A GEANT simulation proved to underestimate the 
measured number of charged particles into the chamber by a factor of 2-4. 
A simulation using stand-alone FLUKA (as opposed to the option of us- 
ing FLUKA subroutines inside GEANT) produced a distribution of particles 
much closer to what was actually seen. Further analysis suggested that a 
significant contribution to the hits in the chambers was due to neutrons in- 
teracting in the mylar face of the chamber and knocking out protons which 
were then detected. Such interactions of low-energy neutrons are included in 
the stand-alone FLUKA. We have not been able to obtain the stand-alone 



45 



FLUKA simulation code. 



3.2.2 GEANT trigger rates 

The trigger rates and cross sections for events in the GEANT simulation 
from the minimum-bias PYTHIA input are shown in Table The total 
fraction of events which pass the trigger in GEANT is slightly more than 
50%, and is higher for lead-in runs than for lead-out by about 1% for all 
types of events. Single diffractive events in which the proton is dissociated 
are triggered on about 15% of the time, and about 10% of those with a 
fragmenting anti-proton pass the trigger. For double diffraction, the fraction 
is about 20%. 

Note that if we add the fraction of triggered single diffractive events (26%) 
times the single diffractive cross section from CDF (9.46 mb) to the non-single 
diffractive fraction (79%) times the corresponding CDF cross section (50.87 
mb), the result is that we should trigger on 42.7 mb, which agrees very well 



with the estimate in Sec. 2.2.3 



3.3 DCC generator 

The operational definition of DCC used by the MiniMax collaboration is 
that DCC is a cluster of pions (domain) with neutral fraction / given by 
the l/{2y/f) distribution and with near-identical momenta which is non- 
relativistic in the DCC rest frame. For the studies presented here, the do- 
main size is taken to be on the order of the detector acceptance, and the 
momentum is such that the DCC is aimed at the center of the acceptance 
with a reasonably large pt in the lab frame. 

We assume that the pions are non-relativistic in the cm frame of the 
DCC domain (and therefore refer to the domain as a "snowball"). We take 
the momentum distribution to be Gaussian with mean (p) = and variance 
(p-p) = 3o"p. Since the DCC domain must be large enough to contain 
physical pions, i.e. have dimensions on the order of a few fm, uncertainty 
principle arguments require a small momentum, cXp ~ 50 — 100 MeV. We use 
the value ap = 100 MeV, which gives the pions a relatively large momentum 
in the lab frame, making them easier to detect. If the pions are not too 
relativistic in the cm frame, then the boosted DCC domain is approximately 
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a circular disk of radius in lego space Rdcc ~ crp/vr- For the MiniMax 
detector, with acceptance approximately a circle with radius 0.6, a domain 
with Rdcc ~ 0.7 — 1.0 would be easiest to find, which implies px = 140 MeV. 
This size seems reasonable since the typical radius of a jet is 0.7. 

We also assume that the number of pions in the DCC domain is ap- 
proximately independent of the pseudorapidity of the center of the boosted 
domain. A result of this is that the ratio of the mean energy density of a 
DCC pion to that of a generically-produced pion is approximately constant. 
Then the mean number of pions is given by 



gen 



Pt \ drjdcj) ) 
or, substituting in the geometric values 



gen -. 




.2 

/I n 

gen 



We take {dN / drf) g^n = 6 in the region of the acceptance, pj^"" = 500 MeV, 
and = 1, which gives a mean number of pions (A^tt) = 5.0. For a given 
event, the number of pions is taken from a Poisson distribution about this 
mean. 

Next the charge of the pions is determined according to the l/{2^/f) 
distribution. The neutral fraction in a given event is generated using the 
transformation method, where, if x is a uniform deviate, then / = is 
distributed according to l/{2\/f). A uniform deviate fji is then generated 
for each of the pions (z = 1, Nj,); if yi < /, the pion is defined to be neutral, 
otherwise it is defined to be charged. 

Figures |3.13| - |3.l5| show for the DCC pions: the momentum distributions, 
the distributions in r] and (f) for a domain aimed at r/ = 4.1 and = 0.75, the 
neutral-fraction distribution and the number distribution. Histograms of the 
number of events with given numbers of charged and neutral pions, for both 
total numbers and the numbers of those that enter the acceptance are shown 
in Fig. |3.16|, and can be compared with that from generically-produced pions 



from PYTHIA in Fig. ^Tf[ 

Unless noted otherwise, all further mention of the DCC generator refers 
to the use of the parameters given here. 
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Figure 3.1: Pseudorapidity distribution dN/dr] for charged pions, all charged 
particles, and photons. 
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Figure 3.2: Energy deposited in (a) scintillator and (b) chambers in minimum 
bias events. 
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Figure 3.3: Energy deposited in (a) scintillator and (b) chambers by single 
charged pions. 
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Figure 3.4: Mean number of wire hits per chamber: real (dashed) and 
GEANT with beam supports, etc (sohd). 
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Figure 3.5: Distribution of NHITS: real (dashed) and GEANT with beam 
supports, etc (sohd). 
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Figure 3.6: Pulse height distributions for chambers 1-4, 9-12, real (dashed) 
and GEANT (sohd). 
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Figure 3.7: Pulse height distributions for chambers 21-24: real (dashed) and 
GEANT (soUd). 
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Figure 3.8: Mean number of wire hits in each chamber (sohd) and with the 
beampipe removed (dashed). 
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Figure 3.9: Mean number of wire hits in each chamber (sohd) and with the 
abort pipe removed (dashed). 
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Figure 3.10: Mean number of wire hits in each chamber (soUd) and with the 
chamber frames removed (dashed). 
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Figure 3.11: Mean number of wire hits in each chamber (sohd) and with the 
beampipe, abort pipe, and chamber frames removed (dashed). 
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Figure 3.12: Charged-particle multiplicity dNch/drj for non-diffractive 
PYTHIA (solid) and HERWIG (dashed). 
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Figure 3.13: Momentum (GeV) of DCC pions: Px and Py in lab frame, Pz in 
DCC rest frame, and Pz in lab frame. 
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Figure 3.14: Location of the DCC pions in lego space. 
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Figure 3.15: Distribution of fraction of DCC pions which are neutral 
{Nj^o/Nj^) and total number of DCC pions {N^). 
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Figure 3.16: Number of charged vs number of neutral DCC pions: total 
number produced, and the number which enter the MiniMax acceptance. 
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Figure 3.17: Number distributions of pions produced generically (by 
PYTHIA). 
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OXo 
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All included subprocesses 


72.98 


300000 


155709 
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152069 


0.5069 


f + f ' ^ f + f ' (QCD) 


2.217 


9012 


7494 


0.8316 
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f + f ^ f ' + f ' 
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98 
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88 


0.6769 
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52456 
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51412 
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0.4667 
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g + g ^ g + g 
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78618 


0.9029 


Elastic scattering 


14.75 


60825 





0.0000 





0.0000 


Single diffractive (XB) 


6.194 


25525 


4057 


0.1589 


3815 


0.1495 


Single diffractive (AX) 


6.194 


25367 


2706 


0.1067 


2452 


0.0967 


Double diffractive 


6.839 


28271 


5920 


0.2094 


5336 


0.1887 


Low-pT scattering 


0.1508 


2104 


1280 


0.6084 


1246 


0.5922 



Table 3.1: PYTHIA cross sections and GEANT trigger rates 
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Table 3.2: NHITS of "Michgan chambers" in GEANT (threshold = 20) for 
varying pulse-height and scintillator cuts. 



Chapter 4 

Data Analysis Tools 



4.1 Tracker 

At least three different track-finding programs were written and used by the 



MiniMax collaboration Each has a different algorithm for reconstructing 
tracks from hit wires, but all give similar results. 

The work described here was done using the combinatorial tracker. That 
algorithm is constructed to find all possible combinations of the hit wires in 
four "crosshair chambers". At least three of these chambers have different 
orientations, so that a unique straight line can be drawn through any such 
combination of hits. The line is considered a potential track, and the non- 
crosshair chambers are searched for hit wires within 3.5 wire spacings of the 
line. If enough wire hits are found, a straight line is fit to the hits. Then, if 
the fitted track passes quality cuts such as cuts on the of the fit, the track 
is recorded. Several sets of crosshair chambers are used in order to increase 
the probability that a track and all the associated hit wires are found. 

Tracks with a greater number of wires are always looked for first because 
they are more likely to be the correct track than a similar one with fewer 
wires; when two tracks are similar (have a large fraction of wires in common) 
the track found first has more weight in determining the track parameters. 
Also, the tracker code is written to search for tracks which are pointed such 
that they appear to originate near the collision point before allowing the 
tracks to point in any direction. Tracks are searched for in the following 
order. First the charged tracks that go straight through the detector are 
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found and recorded. Then the tracking is run on the back 16 chambers and 
tracks which do not share too many wires (the exact numbers are given in 
Table [4.1[ ) with previously-found charged tracks are recorded. These tracks 
include photon conversion tracks and segments of charged tracks which bend 
in the lead and are therefore not found as through-going charged. Tracks of 
the same type are then looked for in the eight chambers directly behind the 
lead (chambers 9-16) in case the tracks leave the acceptance before hitting 
enough chambers. Interactions in the beam pipe also produce background 
tracks which can be found in this region. Finally, the front eight chambers are 
searched for segments of charged tracks that are not found as going straight 
through all chambers. Figures [4.1[ - [4.6| show event displays of various types 
of tracks. 

In the process of searching for each type of track, track candidates are 
compared to previously-found tracks to make sure that the same real track 
is not being recorded as many similar tracks. In the charged tracking, for 
example, tracks with hit wires in at least 22 chambers are searched for first. 
After a track is found, the wires hit by the track are stored in a list of used 
wires. As each new track candidate is found, its wires are checked against 
the used wires. If the candidate does not have at least 17 unique wires hit, 
it is dropped. This significantly reduces the number of fake tracks recorded 
by the tracker. The candidate track is also compared to the other tracks 
individually. The hits in each chamber are compared, and if the tracks are 
within two wires of each other in at least 16 chambers, the tracker determines 
that the same track has been found and the hits of the current candidate are 
added to the list of hits in the previously-found track. We refer to this as 
"grouping". Wires from the grouped track are added to the used- wire list. 
Next, tracks which go through only 21 chambers are considered. If such a 
track candidate has more than 17 wires which are in the used-wire list, the 
candidate is dropped, and if it has 16 wires within two wire spacings of hits 
in another track, the two tracks are grouped. Table ET] shows the various 
cuts on number of chambers hit or wires in common for each type of track. 

The tracks are broken into segments in front of and behind the converter 
plane, referred to as "heads" and "tails", respectively. A dst is written, 
which, though traditionally stands for data summary tape, is just a file which 
contains, for each event, the event number, number of heads and number of 
tails, and then for each track segment, the number of wires hit followed by 
the list of wires. The wire number is given by (chamber — 1) x 128 + wire, for 
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chambers 1-24 and wires 1-128 in each chamber. Earher dst formats included 
the NHITS of the event (which can still be retrieved from the run data file) 
and the track parameters determined by a fit (which are now determined 
using separate code, as discussed in the next section). Sample dst entries are 
shown in Figs. [4.8| . The tracker is run once for each actual or simulated 
run, and the dst is used in all further analysis. 

Only events with an NHITS less than 600 have been analyzed by the 
tracker. Higher-NHITS events tend to have large hit densities in the rear 
chambers, which leads to a huge number of potential tails, most of which are 
not real charged or photon-conversion tracks, and the tracker takes a longer 
time to analyze these events. About 6% of the events in lead-in runs have 
NHITS higher than this value. The NHITS is correlated with multiplicity, 
so that the sample of events which were analyzed by the tracker is somewhat 
biased towards lower-multiplicity events. 



4.2 Track fitter 

A separate program is necessary to fit tracks from the wires written in a 
dst since the fit parameters from the combinatorial tracker are calculated in 
the {x, y) coordinate system instead of in (u, v) (this is important because 
the resolution in u is much better than in v in the chambers behind the 
lead) and uncertainties in the track parameters are not reported. The code 
for the track fitter is given in Appendix B. With the track fitter, straight 
lines are fit to the hit wires in a track, parameterized by u{z) = + buZ 
and v{z) = + b^z, and the covariance matrix (which gives the correlated 
uncertainties of parameters [Q) is generated. Uncertainties in the position 



of the track at a chamber are taken to be 1/ vT2 [| times the wire spacing for 
all hits. 

Another subroutine (see Appendix B) is used to correct a small problem 
with tracks found by the tracker. Occasionally (in about 10% of the events), 
a tail is recorded which, in a single chamber, includes two wires as part of 
the same track which are separated by several wires; these double hits are 
present in most of the non-u chambers of the track. 

^For a hit to be recorded in a single wire in a chamber, a charged track must pass 
within 1 wire spacing of that wire. The uncertainty in position is then given by (Acc)^ = 
- (x)' = ckx^ - (/o ckx)'' = 1/12, where X is the distance in units of wire spacings. 
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The cause turned out to be that two track candidates have common wires 
in all of the w-chambers, and enough close wires in other chambers for the 
tracker to group the track candidates into one track. In some cases, the 
candidates which are grouped in this way are two real tracks, but more often 
one of the candidates is a type of fake track which will be referred to as 
a "ghost". (A ghost is a track which borrows hits in the u chambers from 
a single real track, and finishes up the track with random hits in the non-w 
chambers such as from pipe-shower secondaries.) The track which is recorded 
is halfway between the two candidates. 

Since the occurrence of these tracks with double hits is fairly rare, espe- 
cially for those made of two real tracks, we do not try to separate out two 
distinct tracks, but rather to find one good track. The code used to correct 
this problem tries all possible combinations of hit wires such that the non-w 
chambers each have only one wire included in the track. Only two of these 
combinations are considered as potential correct tracks: the one which best 
points (in v) towards the collision point, and the one which has the best 
when fit to a line in the v-z plane. Of course, a single combination can satisfy 
both these requirements. 

The pointing in v is somewhat complicated by the fact that tracks go- 
ing through the acceptance are at very small angles from the v = plane. 
Therefore, the uncertainty of the z for which f = is much larger than that 
of z{u = 0). In other words, the uncertainty in the angle between the track 
and the v = plane and that angle are both comparably small. Instead, we 
use the following measure of pointing in v. We find the z for which u = 
for the track, defined as zq. Then, working in the v-z plane, we define 6q 
as the angle between the z-axis and a line drawn from [v = 0, z = zq) to 
the point where the track being considered intersects the lead. The angle 
between the track itself and the z-axis in the v-z plane is defined as 6, and 
the uncertainty in that angle as ag. The measure of pointing is taken to be 
\6o — 6\ /a0, which is small for tracks which point (to within uncertainty) to 
the same z as is pointed to in u. 

Figure shows histograms of for tracks with the best x^, \0o — 6\ / ag 
for tracks with the best pointing, and then these values for good tails (those 
which did not have the double- hit problem). The plot of for good tails 
has an obvious separation of what we believe are real tracks and background, 
which we use to choose the cut < 7. For the pointing in f , we require 
\9o — 9\ /ag < 8. We are most interested in tails which point back to the 
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collision point. (Further discussion of pointing cuts is given in Sec. [4.3| .) 
Therefore, if the track with the best pointing also passes the cut, it is 
kept as a good track. If this is not the case, the track with the best is 
considered, and kept if it passes the pointing cut. In most cases, none of the 
combinations of wires produce a track which passes the pointing cuts, and 
the track is dropped. 



4.3 Vertexer 

Charged tracks from the collision point and photon conversion tracks are 
then reconstructed from the track segments. The vertexer code, given in 
Appendix C, is used to determine the probability that track segments meet 
at a common point at the plane of the lead. A charged track is then defined 
as a head which either "matches" at least one tail or passes cuts discussed 
below, and a photon as any group of tails (or a single tail) not matched to 
anything in front of the lead. In terms of the notation (number of heads 
in match, number of tails in match), a charged track is a (> 0, > 0) and a 
photon is a (0, > 0). The location of the vertex is taken as the mean position 
of the included track segments at the lead, weighted by the uncertainties. In 
order to avoid effects from the edge of the lead, the vertex position is required 
to be within the region defined by 4.25 in < u < 10.25 in, —3 in < f < 3 in. 



and y{u — 7.25 in)^ + f ^ < 4 in, which is roughly the area 1 in from all sides 
of the lead. 

The uncertainty in the position of the tail at the lead is apparently under- 
estimated by propagating the uncertainties in track parameters calculated 
by the fitter. The most likely cause of this is multiple scattering in the 
lead. The mean variance in the u position from this error propagation is 
(tJ^ = 5.5 X 10"'* in^ and in v is cr^^ = 0.030 in^. A new estimate was found 
by histogramming the distance between two conversion tracks from a single 
photon in GEANT at the lead, and taking the standard deviation of the his- 



togram as the uncertainty in position. Figure [4.10| shows these histograms, 
which give the new values = 0.007 in^ and = 0.092 in^. (The order 
of magnitude is really more important than the specific value, since cuts on 
other variables used in the vertexer can compensate for small differences in 
these values.) Calculated uncertainties for a track segment which are less 
than the new values are increased to the new values in the vertexer. 
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Tracks are required to point to within some distance of the mean colhsion 
point in order to remove background such as tracks from beampipe shower, 
and also combinatorial fakes. Figures [4.11| and |4.12| show the z for which 
M = of heads and tails of (1,1) charged tracks, and for single (0,1) photon 
conversion tails. Heads of charged tracks are required to point to the region 
—50 < z < 60, which includes almost all primary charged tracks. Most of 
the tracks outside this region are fakes, although some are decay products of 
neutral particles such as the K^. Tails of charged tracks are not required to 
point so that tracks which multiple scatter in the lead will not be dropped. 
For lead-out runs (Fig. |4.11| ), most (0,l)'s which point to z^SO are tracks 
from interactions in the beampipe, while for lead-in runs (Fig. [4.12|) , many 
are from photon conversions in the lead. In order to cut out the fakes from 
pipe shower, single photon conversion tracks are required to point to —40 < 
z < 50. For photon conversions which produce more than one track, at least 
one track in the vertex must point to 2; < 50. The pointing cut in v is taken 
from the study in the previous section to be \{9o — 9)/ag\ < 8. 

Parameters used in the vertexer code were determined by finding those 
which best reconstruct the PYTHIA and GEANT tracks, which are recorded 
in separate files. The pointing cuts, the increased uncertainties in the u and 
V position of tails, and cuts on the ior the fit of how well two tracks match 
at the lead (which will be discussed in Sec. [4.3.1| , see Fig. [4.13| ) are not 
independent; the cuts chosen are all self consistent. The code seems to work 
sufficiently well for both simulated and real data. As an example. Fig. [4.14 
shows the tracks which survived the vertexer cuts for the large NHITS event 



of Fig. U 



Plots of efficiencies for finding charged tracks and photons are shown in 
Figs. |4.15| - |4.2^ . Efficiencies are given as a function of energy, transverse 



momentum, multiplicity, NHITS, and position in rj and 0. The numbers 
of charged tracks and of photons in the acceptance as a function of these 
parameters are also shown, both because these distributions are interesting, 
and also to give an indication of the statistical significance of bins in the 
efficiency plots. The plots also include the mean number of fakes as a function 
of multiplicity and NHITS, and the mean number of actual charged tracks 
and photons as a function of NHITS. 
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4.3.1 Matching tracks at the lead 

Tracks are "matched" or "vertexed" in the following way. Consider two tracks 
intersecting the lead at points = = {ui,Vi), ^2 = (^2^, ^2^) = 

{u2,V2), with covariance matrices Ci and C2. Assume that the two tracks 
meet at exactly the same point [$, = {^^,^'^)] at the lead. Also assume that 
the are Gaussian distributed about mean ^. 
The for the fit of the two tracks to ^ is 

x' = (li - i)Ci-\^i - + (6 - ^)C2-^(6 - 1)- (4.1) 

Minimizing with respect to ^ * gives 

C = Rcr' + [(cr'y'^i' + {02-^%"] . (4.2) 



Putting this ^ back into Eq. [4.1| yields x^- There are four known parame- 
ters (^1 and ^2), and two which are determined (^), leaving two degrees of 
freedom. The reduced is therefore 

All track segments in an event were matched together in order to find the 
cutoff in for which matches with a lower are most likely real tracks, 
while those with a higher are most likely unrelated track segments. Plots 
of matches between pairs of heads, pairs of tails, and head-tail pairs are shown 
in Fig. |4.13| . Head-head matches which have a low x^ are due only to the 



heads being coincident ally close together at the lead. The peak at low x^/2 
for head-tail matches is due to charged tracks, and is fairly cleanly separated 
from false matches. The tail-tail plot does not show a clear division, but a 
slight division is present at about the same value of x^/2 as the division for 
head-tail matches. (Note that the cutoff values could have been different for 
the head-tail and tail-tail matches since the uncertainties in position of the 
tails are much larger than those for heads.) A reasonable cut appears to be 
log(xV2)<0.7, or xV2<5. 



4.3.2 Tracks which may be grouped together 

The original vertexer algorithm grouped any track segments which matched 
any other track segments in a vertex. The uncertainty in position of the tails 
is large enough that this occasionally produced vertices which contained more 
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than one real track. For instance, a (1,1) charged track and a (0,n) photon 
could be grouped into a (l,n+l), or two charged tracks could become a (2,2). 

By demanding that all tails in a vertex match the head (which has a much 
smaller uncertainty), charged tracks are less likely to be grouped with other 
tracks. The vertexer code was therefore changed to look for charged tracks 
first and remove those tails which were matched to heads before matching 
tails together. The following is the result of this change for a sample of 
GEANT events. 

Of the original 135 (2,2) vertices, 104 have heads which do not match 
each other and were therefore separated. About 94% of these are really 
two charged tracks, and 6% are fake heads grouped with real tracks. The 
31 vertices which have two heads that match each other consist of 48% fake 
heads along with other tracks, 13% e''"e~ pairs from photon conversions in the 
window ("window conversions"), and 39% vertices with two charged tracks 
which are very close together at the lead. 

The majority of (l,n)'s which were separated into more than one vertex 
are charged tracks combined with tails from photons, fake tracks, or pipe 
shower tracks. The latter two types of tails are likely to be removed by 
pointing cuts when separated from the charged track. About 10% result 
from window conversions where a conversion secondary showers in the lead, 
and approximately 5-10% are single charged tracks with multiple tails (from 
interactions in the lead). Of those which remained (l,n)'s, approximately 
20% are pions which interact in the lead, 30% are pairs of a charged track 
and a fake tail, 25% result from window conversions, and 20% are pairs of a 
close charged track and a photon. The rest involve missing heads, decays in 
GEANT, etc. Only about half of the vertices which contained both a charged 
track and a photon were reclassified as such. 

The overall improvement in vertexer performance due to these changes 
is small since these types of vertices are not very common. Efficiencies are 
improved for high multiplicity events, without introducing many fakes. 

Studies were also done with the real data by adding a single charged track 
or photon conversion track from a clean event to every other event in the run. 
This was done at the level of the dst, so that the behavior of the tracker for 
the combined events is ignored. At this level, the efficiencies for finding all 
tracks of the individual events in the combined event are not compromised 
unless the added track is within some distance of a track in the other event. 
The mean distance between an added charged track and another charged 
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track for a track to be missed is about 0.5 in, and for an added charged track 
and a photon conversion about 0.7 in. 

4.3.3 Tail-less heads as charged tracks 

Often the tail of a low-energy charged track will not be detected due to 
multiple scattering in the lead which may cause it to stop or to bend out of 
the acceptance or to lose so much energy that it does not travel in a straight 
path through the rear chambers. In GEANT, most of the (l,0)'s are real 
(collision-point) charged tracks which stopped in the lead or had very soft 
tails; some are from neutral particle decays. Those tail-less heads which are 
fakes are distinguished by the fact that they do not have hits in all eight of 
the front chambers. However, in the real data, the (1,0) 's seem to include a 
large fraction of fakes. The ratio of the number of (1,0) 's to the total number 
of identified charged tracks [including (l,0)'s] is much higher for the real data 
(run 1125) than for GEANT, 17% vs 7%, whereas if the energy distributions 
of charged tracks are similar for real data and GEANT, this fraction should 
be about the same. Ehminating those (1,0) 's which do not have hits from 
all eight chambers reduces the fractions to 12% and 6%. Other types of 
(l,0)'s which appear to be fakes either share wires in all three u or all three 
V chambers with another track (ghosts), or point to somewhere other than 
the intersection with the 2;-axis (in u) of other charged tracks in the event 
which go all the way through the detector. 

Based on these observations, tight cuts are made which are intended to 
eliminate most fakes at the expense of losing some real charged (l,0)'s. Any 
(l,0)'s which share wires in all u or all v chambers arc dropped. The pointing 
cut is defined by other charged tracks in the event; if there are none, it is not 
used. The collision point Zcp is determined as the weighted mean of the point 
of intersection with the 2;-axis of any other heads of (1,> 0) charged tracks 
in the event, using the parameterization in u of the heads and the associated 
uncertainties. Plots of the distance between the z{u = 0) of the (1,0) and 
the mean Zcp divided by a, the root mean square of the uncertainty in these 
z, for both GEANT and real-data (l,0)'s are used to choose the condition 
that the z{u = 0) of the (1,0) must be within 2a of Zcp- This cuts out many 
real charged (l,0)'s. These cuts reduce the ratio of "charged" (l,0)'s to total 
charged tracks to 5.0% for GEANT and 9.6% for run 1125. Although the 
fraction for the data is almost twice as large as that for the GEANT, no 
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apparent qualities of the remaining (1,0) 's suggest that they are fake. 

4.3.4 Middle-eight tracks 

Originally, tails which do not go through at least 14 of the rear chambers, 
but do go through the first eight chambers behind the lead were not used by 
the vertexer. The reason is that the rear chambers are often flooded with 
pipe shower, which results in a large number of "middle-eight" tracks which 
arc pipe shower tracks or combinatorial fakes which can be found as photon 
tracks by the vertexer. However, these tracks can be used carefully in certain 
circumstances. Allowing the middle-eight tracks to be vertexed as tails of 
charged tracks saves some real charged tracks that would otherwise be (l,0)'s 
and might have be thrown away by the tight cuts. The vertexer also allows 
middle-eight tracks as photon conversions as long as they are vertexed with 
at least one other tail. This is important because the tight pointing cuts on 
(0,1) photons are often failed by a single tail of a photon conversion which 
produces other tracks that leave the acceptance before reaching the final 
chambers. If a middle-eight tail is vertexed to such a tail then the vertex is 
not subject to the tight cuts. Also, since fake tracks are unlikely to vertex at 
the lead, a vertex made only of at least two middle-eight tracks is counted 
photon. 

It turns out that very few events have middle-eight tracks which can 
be vertexed with other track segments at the lead, so that the effect on 
overall efficiencies is almost negligible. This may be related to the difficulty of 
getting a good fit for middle-eight tracks due to the lack of f -resolution. Fake 
charged tracks are occasionally created when a fake (1,0) is coincidcntally 
matched with a middle-eight tail. The middle-eight tracks were nevertheless 
included in the vertexing, since even a tiny improvement in the photon- 
finding efficiencies is welcome. 

4.3.5 Origin of fakes 

A charged track from the vertexer is classified as a fake if there are no charged 
tracks from PYTHIA inside the acceptance within the specified distances in 
u and V of 0.5 in and 1 in, respectively, of the identified track at the lead. A 
sample of events which have one fake charged track was taken from 1.5 x 10^ 
GEANT events. In order to avoid effects due to the presence of photons. 
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the events were required to have no photons observed or known to convert in 
the acceptance. The sample is divided into the following categories: events 
in which 0, 1, or 2 charged tracks are sent into the acceptance and one 
extra charged track is found (the total sample of events had 3402, 1475, 348 
events, respectively, in these classes, and 39, 38, and 38 events were used in 
this study), and those in which charged tracks are sent into the acceptance 
and 2 are found (217 total events, 37 events used, accounting for 74 fake 
tracks). Of the 188 fake charged tracks, 
150 are decay product (s) of K^'s or A's, 

24 are secondary charged tracks from other decays or interactions in the 
material surrounding the detector 

1 is from the conversion of a photon in the window, 

12 are real charged tracks just outside the acceptance, and 

1 is an actual combinatorial fake. 
Decays of single K^'s or A's are responsible for about 70% of the events with 
two fake charged tracks. The important outcome of this study is that almost 
all of these "fakes" are caused by interesting physics processes or edge effects. 

Also studied were 20 events (from a total of 416) where a photon conver- 
sion is present in GEANTQ and that photon is not found, and no PYTHIA 
charged tracks are sent into the acceptance, but one is found. In 11 of these 
events, the photon converts in the window, and a resulting e^ showers in 
the lead, i.e., the missed photon and fake charged track are correlated. The 
remaining events involve photons which are not found by the vertexer or 
are found just outside the acceptance, with uncorrelated fake charged tracks 
which fall into the categories mentioned above. 

A photon found by the vertexer is defined as fake if no PYTHIA photon 
is aimed such that it will hit the lead within 0.5 in in u and 1 in in f of the 
identified conversion. In 38 events (from a total of 1800) where no charged 
tracks are sent into the acceptance by PYTHIA and no photon conversions 
from GEANT are present in the acceptance, and one photon is found, 

19 are conversions of secondary photons, roughly half of which appear 
to be due to decays, and the other half due to interactions in the 
detector region, 

^ GEANT conversions are defined as e^ tracks originating in the region containing the 
lead and scintillator between chambers 8 and 9, and within 0.1 in at the lead of a photon 
from PYTHIA. 
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14 are real photons just outside the acceptance, 

2 are charged decay products of K^'s where the heads are lost, 

1 is the product of an interaction of a neutron in the lead, and 

2 are from events which are too complicated to determine the source 

of the fake. 

Of 20 events (from a total of 319) where a charged particle is lost and a 
fake photon found, 17 involve a charged track which multiple scatters in the 
lead so that the tail is found as a photon conversion, but the head is dropped. 
In the remaining 3 events, the missed and fake tracks are uncorrelated. 

When a pipe shower produces a large density of hits in the rear chambers, 
the probability is greatly increased for finding fake photons. Combinations 
of these hits, possibly together with hits from real tracks from the collision, 
form tracks which point to the collision region and therefore should be found 
by the tracker. Ghosts (described in Sec. [4.2| ) are a common example of such 
a fake track. 

If pipe showers were uncorrelated with important observables such as the 
number of charged tracks and photons entering the acceptance in an event, 
then vetoing events which have pipe shower would reduce the number of fake 
photons found without biasing the remaining sample of events, say towards 
those with low multiplicity. Unfortunately, this is not the case. In order 
to study this, we first identified events with pipe-shower tracks. Primary 
tracks which are produced in the collision and pass all the way through the 
MWPC telescope occupy a very limited region of the phase space defined 
by the position and angle in u at some z. (Recall that the acceptance is 
roughly defined by the coverage of the lead converter, which extends from 
u = 3.25 in to 11.25 in and v = —4.0 in to 4.0 in at ziead = 150 in.) Figure 



4.25| shows the location of tracks in {u, 9u) space where u = u{ziead) and 
tanOu = u{ziead) / ziead, averagcd over many events. The dense band near 
the center of the plot is due to primary charged tracks. We defined rough 
boundaries around this region and vetoed events with any tails in the regions 
9u > 0.009U + 0.016 or 9u < 0.006u- 0.028, excluding tails which the vertexer 
found as photon conversion tracks. (Note that this is a very tight cut.) 

Because of the higher hit density in the rear chambers created by pipe 
shower, the events which remain after the pipe-shower veto have a much 
lower mean NHITS than that for all events. In GEANT, the mean NHITS 
for lead-in runs dropped from 142 to 90, and in run 1125 (lead in) from 229 to 
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113. Also, surviving events with higher NHITS tend to have more real tracks, 
as is shown for GEANT in Fig. [4.26 . This figure also shows the distributions 
of actual charged-track and photon multiplicities (from PYTHIA), which 
appear to be biased towards lower multiplicities. This seems to be true to 
a greater extent for the real data; Figure [4.27| shows that the probability 
distribution of the observed number of photons before and after the pipe- 
shower veto for real data is more biased towards lower multiplicities than 
that for GEANT. The mean number of charged tracks found decreased from 
(^c/i) = 0.53 for all GEANT events to 0.47 for this sample, and of photons 
from {n^) = 0.21 to 0.14. For the data (run 1125), the cut is even harder: 
{nch) = 0.48 falls to 0.32 and (n-y) = 0.19 drops to 0.06, practically no 
photons. We expect the veto to have a greater effect on the data because 
we believe that the real events contain more pipe shower than the GEANT, 
but since the vetoed events are apparently correlated with high multiplicity 
events, the remaining events are a clean, but not unbiased, sample. 

This sample of non-pipe-shower events is useful for studying fake photons 
which were not caused by pipe shower, although the events are also less likely 
to contain converted photons. Fake photons in two sets of the remaining 
events were classified. The first set was taken from about 10^ minimum-bias 
(triggered) events, and contains 104 events which pass the pipe-shower veto 
and have a fake photon. The second contains 71 such events which also have 
at least four primary charged particles and/or photons (from PYTHIA) in the 
acceptance; this sample includes all events of that type in the approximately 
1.5 X 10^ triggered events. 

Of the fakes in the first set, 

5 are ghosts, 

11 are tails of charged tracks which are not vertexed with the heads, 
9 are associated with photon conversions but are not vertexed with 
the other conversion tracks, 

31 are other secondary photons such as those produced by interactions 
in the beampipe, 

2 are secondaries from the interactions of neutrons in the lead, 
11 are photons from the decays of 7r°'s produced by — > 7r°7r°, 
18 are real photons just outside the acceptance, 
and the remaining 17 either are not able to be classified, or fell into unique 
categories. The first few types of fakes are the most serious because they 
are correlated with the presence of real tracks. In almost all cases, when the 
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tail of a charged track is not vertexed to the head, the tail-less head does 
not pass the cuts necessary to be counted as a charged track. Therefore a 
charged track is lost and is replaced by a fake photon. (This comphcates the 
generating function analysis described in Chapter 5, because it introduces 
a correlation between charged-track and photon efficiencies.) The classifica- 
tion of photons which are outside the acceptance but arc found inside the 
acceptance as fakes is really just due to details in the code for finding fakes. 

Due to the higher multiplicities in the second set, we expect more fakes 
associated with real tracks, such as ghosts. We find that this is true, but not 
to an extent which would seriously invalidate the assumptions made in the 
next chapter. The fake photons include 

8 ghosts, 

25 tails of charged tracks, 

12 tracks associated with a photon conversion, 

11 other secondary photons, 

2 photons from n^n^, 

9 photons just outside the acceptance, 
and 4 unclassified fakes. 
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Figure 4.1: Event display of a clean charged track. 
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Figure 4.2: Event display of a straight-through charged track and a char 
track which bends in the lead. 
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Figure 4.3: Event display of a photon with three conversion tracks. 
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Figure 4.4: Event display of two clean photon conversions. 
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Figure 4.5: Display of an event with two charged tracks and a photon, and 
typical NHITS. 
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Figure 4.6: Display of an event with a large NHITS. 
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Figure 4.7: Entries in the dst for the events shown in Figs. 4.1 and 4.3. 
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Figure 4.8: Entry in the dst for the event shown in Fig. 4.5. 
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Figure 4.9: Histograms of for the fit in v and of \9q — Q\ joQ for re-fitted 
tails and good tails. Cuts were made of < 7 and \Qq — Q\ j gq <'^. 
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Figure 4.10: Distance in u and v between photon conversion tracks at the 
lead when there are two conversion tracks in GEANT. 
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Figure 4.11: Pointing, defined as tlie z wliere the track goes through w = 0, 
for heads and tails of charged tracks, single photon conversion tracks, and 
the z where any two heads have the same u position, for lead-out run 1093. 



91 




Figure 4.12: Pointing, defined as tlie z wliere the track goes through m = 0, 
for heads and tails of charged tracks, single photon conversion tracks, and 
the z where any two heads have the same u position, for lead-in run 1092. 
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Figure 4.13: Values of log \ for head-tail, head-head, and tail-tail matches at 

2 2 

the lead. A superposition of these plots shows that log ^ = 0.7, or ^ = 5 is 
a reasonable cutoff between matching real tracks together and track segments 
which are coincidentally close at the lead. 



93 





/ 

/ 

1/ 


1 


y • . • : X 




V . i X 

/ % / 














A'- i 






* ! * • i ^.^--^^^^^-^^^^ 
? . • : » » : : 




: • . - . ! - ; 











O 25 50 75 100 125 150 175 200 225 250 275 



evnum= 1363, nhits= 578 



Figure 4.14: Display of the event shown in Fig. 4.6 for tracks kept by the 
vertexer. 
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Figure 4.15: Efficiencies for finding charged tracks as a function of energy 
and transverse momentum of the track, number of charged tracks vs energy, 
transverse momentum. 
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Figure 4.16: Efficiency for finding piiotons as a function of piioton energy, 
efficiency for finding photons know to convert vs energy, probabifity of con- 
version of a photon vs energy, number of photons vs energy. 
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Figure 4.17: Efficiency for finding pliotons as a function of photon transverse 
momentum, efficiency for finding pfiotons know to convert vs pr, probability 
of conversion of a photon vs pr, number of photons vs pr- 
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Figure 4.18: Efficiency for finding cliarged tracks as a function of total multi- 
plicity into the acceptance, efficiency for finding photons vs total multiplicity, 
frequency of observing an event with a given total multiplicity. 
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Figure 4.19: Efficiency for finding cliarged tracks as a function of charged 
multiplicity into the acceptance, efficiency for finding photons vs photon mul- 
tiplicity, frequency of observing an event with a given charged multiplicity, 
frequency vs photon multiplicity. 
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Figure 4.20: Mean number of fake charged tracks found per event with a 

given total muhipUcity, mean number of fake photons found per event with 
a given photon muhiphcity, mean number of fake charged tracks vs charged 
multiphcity, mean number of fake photons vs photon multiphcity. 
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Figure 4.21: Efficiency for finding charged tracks as a function of NHITS 
of the event, efficiencies for finding converted photons and all photons vs 
NHITS, frequency of observing an event with a given NHITS. 
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Figure 4.22: Mean number of fake charged tracks found per event with a 
given NHITS, mean number of fake photons found per event with a given 
NHITS, mean numbers of real charged tracks, real photons given NHITS. 
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Figure 4.23: Efficiencies for finding cliarged tracks, photons at a given rj, 
number of cliarged tracks, photons vs rj. 
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Figure 4.24: Efficiencies for finding charged tracks, photons at a given (j), 
number of charged tracks, photons vs 4>. 
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Figure 4.25: Location of tracks at the lead in phase space {u, 9^) for lead-in 
GEANT, lead-out GEANT, lead-in run 1110, lead-out run 1111. 




Figure 4.26: PYTHIA charged-track and photon multiplicity distributions, 
mean multiplicities vs NHITS for all lead-in GEANT events (solid) and for 
non- vetoed events (dashed). 
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for non-vetoed events from lead-in run 1125 (solid) and lead-in GEANT 
(dashed) . 
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Table 4.1: Cuts used by the tracker. 



Chapter 5 



Generating Function Formalism 
and Robust Observables 

In order to determine the distribution of the neutral fraction (/) of pions in 
the MiniMax detector, we would like to count the number of charged and 
neutral pions from a collision entering our acceptance. (Even if only a piece 
of a DCC domain enters the acceptance, we would expect on average to 
observe the same / as for the entire domain.) This approach is complicated 
by many things, including the fact that vr^'s decay almost immediately into 
two photons which we identify only if they convert in the lead or scintillator 
within the MWPC telescope. To precisely determine the number of vr'^'s, we 
would need to identify both 7's from a decay and reconstruct the pion mass 
using the electromagnetic calorimeter. 

We are unable to do this due to the fact that the probability of both pho- 
tons from a decay entering the small acceptance of the calorimeter is only 
about 15%. However, we have found a set of observables which sample the 
charged-neutral distribution and are independent of many of the detector- 
related complications and (uncorrelated) efficiencies, and take very different 
values for pure DCC and for generic particle production. We make many 
bold assumptions about the production and detection of particles going into 
the detector in order to establish the robustness of these observables; how- 
ever, simulations indicate that the assumptions are not unreasonable (see 
Sec. |6.1.1|) . We assume that particles other than pions can be ignored, that 
charged particles and photons are not misidentified, that the production pro- 
cess can be modeled as a two-step process, with the total number of pions 
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given by a parent distribution and the fraction of pions that are neutral given 
by, e.g., a binomial or DCC distribution function, and that detection efficien- 
cies for finding a charged particle or photon do not depend on the nature of 
the rest of the event. 

These observables and their properties are best understood using gener- 
ating functions and their factorial moments to describe probability distribu- 
tions for the production of some species of particle. The generating function 
formalism has been widely used in multiparticle analysis pT|-p5| and was ex- 



tended to two variables |^-||5^, [@, e.g., charged and neutral pions by the 
MiniMax Collaboration in Ref. |^. The ideas in this section are described 
in detail in that paper. 

Some of the difficulties in measuring the charged-neutral distribution are 
listed below: 

1. The MiniMax acceptance is small, so that it is improbable that both 
7's from a 7r° enter the detector acceptance; 



2. the conversion efficiency per 7 is only about ou7o.; 

3. not all 7's come from vr^'s; 

4. not all charged tracks come from vr^'s; 

5. because of the small acceptance, the multiplicities are rather low, so 
that statistical fluctuations are very important; 

6. detection efficiencies are not the same for charged tracks and 7's and 
are momentum-dependent; 

7. efficiency functions may depend on the observed multiplicity or other 
parameters; 

8. the efficiency for triggering when no charged track or converted 7 is 
produced within the acceptance is relatively low and different from 
that for events in which at least one charged particle or converted 7 is 
detected. 
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5.1 Generating functions for charged-pion— 
neutral-pion distributions 

If the set of (normalized) probabilities for producing N particles in a given 
region of phase space is {P{N)}, then the generating function can be defined 
as 

oo 

Giz) = E ^""m), (5.1) 

N=0 

and contains all the information of the {P{N)}: 

Information can likewise be extracted from the factorial moments, defined as 

f^^[^^) ^ = {N{N-l)...{N-z + l)), (5.3) 

where (C) = EnOP{N). 

Now if p(ncfe, no) is the probability distribution for producing rich charged 
and no neutral pions, the generating function is 



oo oo 



G{zch:Zo)^ J2 T.Pi^ch:no)z:tz^\ (5.4) 

"ch,=0 no=0 

and the factorial moments for charged (ch) and neutral (0) pions are 

We assume that p{nch,n'o) can be written as the product of a parent 
distribution P{N) for producing total pions, and p(nc/i, no; N), which gives 
the charged- neutral distribution of pions for a given N — rich + ^o- 

P{nch: no) = P{N)p{nch, no; iV), (5.6) 

where 

oo 

Ar=o 



Ill 



oo oo 



J2 ^N,ncn+no pinch, no; N) = 1. (5.8) 

nch=0 no=0 

For generic production of pions, the charge is distributed according to a 
binomial {bin) with mean neutral fraction / = 1/3, so that 

PMn{nck,no;N) = -—^ -/""(I - (5.9) 

no!(iV - no)! 

The corresponding generating function [from Eqs. ( |5.9| ), ( p. 61) and (|5.4| )] is 

GMn{zch,zo;f) = Y.PiN)[ho + {l-f)z,hf. (5.10) 

N 

Note that Gfej„ depends only on the linear combination 

C = fzo + (1 - />ch. (5.11) 

Much of the simplicity of the generic case is also realized for what can be 
called the binomial transform, 

p{n,H,no;N) = -—^ [' dfp{f)r{l - fr^\ (5.12) 

no!(A' — no)] Jo 

of the arbitrary normalized distribution p{f), such as the DCC distribution 
p{f) = 1/(2 v7). This leads to a wide class of possible pion factorial-moment 
generating functions, namely 

G{Zch,Zo)= dfp{f)Gun{ Zch, Zo', /), (5.13) 
Jo 

where Gbin{zch, zo] f) is given by ( [5.10| ) with / replaced by an arbitrary /, 
0</<l. 

If P{N) is a Poisson distribution, [P{N) = ^e"^ with /i = (A^)] then 

Gbin{zch, Zo; f) = e~^+''^ , so that In Gbin{zch, zo] f) is linear in (. The PYTHIA 
simulations yield generating functions that, to good approximation, depend 
only on a fixed linear combination of Zch and zo (Fig. |5.1j ); the full detector 
simulation with GEANT is found to alter this linear behavior slightly. Com- 
pare this to the generating function for the DCC distribution which depends 
on both Zch and Zo (Fig. |^ ). 
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5.2 Generating functions for charged-pion— 
photon distributions 

Next we take into account the probability ech for observing a given primary 
charged pion in the detector and a probabihty (1 — ech) ioT not observing 
it, and the probabihties e^, m = 0, 1, 2, for observing m photons from a 7r° 
decay, with eo + ei + e2 = 1. We assume that these efficiencies are uncorre- 
lated. Then the generating function for the distribution of observed particles, 
including efficiencies, is obtained from G{zch, zq) [0 ^y replacing Zch by the 
generating function 

gch{Zch) = (1 - (^ch) + ^chZch, (5.14) 

and Zq by the generating function 

g^iz^) = eo + eiz^ + €2^7^ (5.15) 

For the class of production models characterized by ( |5.13| ), this leads 
to the following factorial-moment generating function for the distribution of 
observed charged pions and photons: 

Zch ) ^7 / / dfp{f)Gbin{gch{zch),go{z^)', f)- (5.16) 

J 

The charged-pion-photon factorial moments are 

which introduces the indexing (^, j) with respect to charged particles and 
photons which will be used in the remainder of this section. The two lowest 
orders of factorial moments are 

/i,o= M ={l-f)ech{N), 
hi= K) =(/)(ei + 2e2)(iV), 

/2,o = {richinch - 1)) = ((1 - f?) ech' {N{N - 1)), (5.18) 

= (richn,) = (/(I - /)) echiei + 262) {N{N - 1)), 

/o,2= {n,{n,-l)) =(/2)(ei + 2e2)2(iV(iV-l)) + (/)2e2(iV), 

where the overall statistical averages for the charged, photon, and charged- 
photon factorial moments are expressed, in an obvious notation, in terms 
of the independent moments taken with respect to the P{N) and p{f) dis- 
tributions. The second-order factorial moments represent the lowest-order 
correlative effects among charged pions and photons. 



113 



5.3 Robust observables 

We would like to construct a measure from the moments in the form of a 
ratio in order to cancel out as many effects as possible, apart from the p{f) 
averages which give information about the charged-neutral distribution. The 
gamma-gamma correlation (/o,2) involves a term proportional to €2 (N) that 
cannot be cancelled by any other moments of this order. However, the ratio 

^ ^ /i,i/i,o ^ (nchn-y) {rich) . . 

^'^ /2,o/o,i {nchiuch - I)) (n^) ' 

involving the other four moments has complete cancellation of all reference 
to the background distribution P{N) and the efficiencies ei, €2, and ech for 
generating functions of the form ( ^.16| ): 

,^ (/(I -/))((!-/)> 

- {(i-m(/> ■ 

For generic pion production, p{f) = 6{f — /) and 

n^iigen) = 1. (5.21) 
For a DCC distribution, p{f) = l/{2^/f), the ratio is 

n^iiDCC) = ^ . (5.22) 

Therefore, the pure DCC and generic distributions should be easily distin- 
guishable if the statistical uncertainties are not too large. 

We can simplify the formulas slightly by introducing normalized factorial 
moments. For the production of particles, these are given by 

A generalization of the Fj's to normalized moments for charged track and 
photon production is 

_ {richjnch - 1) ■ ■ ■ {uch - ^ + 1) %(^7 - 1) ... (n^ - J + 1)) 

{rich) {n-r) 
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In particular, 

_ F, ((1 - /)•> 
•'° " {(I-/))' * ' 

and 

^ ^ -P^m (/(I - /)') 

(/)((! -/)r ' ^ ^ 



where Fi refers to the ith normahzed factorial moment ( |5.23| ) of the P{N) 
distribution for the total multiplicity. 

We see that ri^i = Fi^i/ F2fi. A generalization of ri^i to a family of robust 
observables is 

, _ ^.1 _ ((1 -/))(/(!-/)-) ..27) 

where again the dependence on the parent distribution and efficiencies has 
dropped out. For alH > 1, generic particle production yields 

r,^i{gen) = 1, (5.28) 

while for DCC, 

n^i{DCC) = . (5.29) 

% + 1 

Thus, Tj 1 becomes more sensitive to the difference between DCC and generic 
production mechanisms with increasing order of the moments. This reflects 
the broadness of the DCC neutral-fraction distribution relative to the bino- 
mial distribution of the generic case. If only a fraction of particle production 
is due to DCC, the signal will be easier to see in the higher-order ratios, 
which are sensitive to the tail of the charged-photon distribution, where the 
ratio of DCC to generic production is relatively high. The ratios 

= -p^ (5.30) 

involving higher-order gamma correlations (-Fij, j > 1) are not robust be- 
cause the moments Fi j depend on the photon detection efficiencies. However, 
the moments can be expressed in terms of only one combination of these ef- 
ficiencies along with the mean number of photons, namely 



(5.31) 
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as 

[i/2] 



^ij — 5Z ^i/'"^ -^«+i-ni//T £\\i I £\j-m ' (5.32) 



The coefficients Cj^m are obtained from the identity true for any differ 
entiable function Diz^), 



^ ""'^ ' E c,,^2-{2zy-^-^^J^ . (5.33) 

m=0 



The first few Cj,m are 

Cj,0 = 1, 

c,,i = j(j-l)/2, 

c,,2 = 3 j!/4!(j-4)!. (5.34) 

The ratios can be used in the analysis of experimental distributions, with 
the understanding that the parameter ^ is to be determined from the data. 
Generally, we have the bounds and limiting values 

rij{gen) > 1, (5.35) 

ri,ji9en)\^^Q = 1, (5.36) 



^!(2j - 1)!! 

Finally, we turn to the effect of the MiniMax trigger on the moments 



nADCC)\^=,= ■\^\ J- . (5.37) 



and ratios. As discussed in Sec. |2.2| , events in which no charged particle or 
converted 7 goes through the acceptance of the detector are triggered with 
lower efficiency, e, than other events. Take p"^^{nch, n.y) to be the (normalized) 
probability for observing an event with rich charged particles and converted 
7's in the acceptance assuming perfect triggering, and p*'^^^{nch, n^) to be the 
probability including the effects of both the trigger and the particle detection 
efficiencies. An effective model for the effect of the trigger on the probability 
is given by 

p*"9(0, 0) = e«p°'^(0, 0), Tick = = 0, (5.38) 

p^"^{nch, n^) = ap°^\nch, n^), n^h + > 0, (5.39) 
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where a is a normalization factor, 



1 + (1 - e) p°''"(0, 0) 



(5.40) 



When we incorporate these trigger efficiencies, 




(5.41) 



and 




a 




(5.42) 



leaving the r^^j robust in the sense that they are also independent of e. 
In summary, 

1. the do not depend upon the form of the parent pion multiplicity 
distribution; 

2. the J arc independent of the detection efficiencies for finding charged 
tracks, provided these efficiencies are not correlated with each other or 
with other variables such as total multiplicity or background level; 

3. the Ti^i are also independent of the 7 efficiencies in the same sense as 
above; the ?"j,(j>i) depend only upon one parameter, S,-, which reflects 
the relative probability of both photons from a 7r° being detected in 
the same event; 

4. in all cases the r^^j are independent of the magnitude of the null trigger 
efficiency; 

5. the ratios rj j possess definite and very different values for pure generic 
and pure DCC pion production. 
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Figure 5.1: Contour plots of the generating function G{zch,Zo) and 
log G (2:^/1, -2o) from all PYTHIA charged and neutral pions, and from those 
entering the MiniMax acceptance (the vertical axis is Zch and the horizontal 
is Zq). 
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Figure 5.2: Contour plots of the generating function G{zch, Zq) and 
logG{Zf^^^fij Zq) from all charged and neutral pions from the DCC generator, 
and from those entering the MiniMax acceptance (the vertical axis is Zch and 
the horizontal is 2:0). 



Chapter 6 
Results 



6.1 DCC: Calculated factorial moments and 
robust observables 

For each real or simulated event, the vertexer returns the number of charged 
tracks and photons found. Then the frequency distribution for observing a 
given number of charged and photon-conversion tracks, J\f{nch, n,^), is used to 
calculate the factorial moments, assuming that the probability of observing 
an event with a given number of charged tracks and converted photons is 
given by the observed distribution: 

P{nch, n^) ~ ^f{nch, n^)/N, (6.1) 
where N is the total number of events 

^= E E^f{nch,n,). (6.2) 

nch=0 n-f=0 

Statistical errors are estimated assuming Poisson fluctuations and stan- 
dard propagation of errors formalism. Appendix A shows the derivation of 
some of the formulas used to calculate uncertainties in the moments and rij. 
In order to check the accuracy of the calculated uncertainty in the rjj, the 
values of ru and 0",.^ were determined for groups of 25000 events from nms 
1096, 1099, 1103, 1109, 1110, 1125, 1126, 1127, 1137, and 1139. Then a 
Gaussian was fit to the histogram of (r(i) — {r))/ar{i), where r{i) and ar{i) 
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are ri^i and <Jri,i, respectively, for the ith group of events, and (r) = 1.0228 



is the mean ri^i for all events. The fitted Gaussian is shown in Fig. 3.1 



and has mean /i = 0.012 ± 0.144 and standard deviation a = 0.985 ± 0.160. 
Note that if the calculated uncertainty is equivalent to the actual statistical 
uncertainty, then the standard deviation of the Gaussian should be cr = 1. 
Since this is true to within errors, we believe that the calculated uncertainty 
accurately represents the true statistical uncertainty. 



6.1.1 Simulations 

The moments and rij were calculated for approximately 1.5 x 10^ PYTHIA 
events which would be seen by the detector (pass trigger cuts) and 2 x 10^ 
pure DCC events, and the results are shown in Table |0|. The PYTHIA 
results are given for perfect charged- and photon-finding efficiencies, along 
with the output of running these events through the GEANT simulation. 
The DCC events are also processed by the same GEANT simulation, except 
that the trigger is not used since there are no particles generated in the 
p direction. For purposes of comparison, the predicted values for ideahzed 
binomial and 1/(2^7) distributions are included. [For j > 1, these depend 
on the parameter ^, which was estimated using the relationship between 
/o,2, /2,o assuming a binomial distribution. Then the observed (n^) from 
PYTHIA/GEANT was used to obtain 2e2/(ei + 2e2) ^ 0.08 ± 0.01. The Fi 
used in these predictions were also determined from the PYTHIA/GEANT 
data. We do not expect these values to be correct for the DCC case; in 
particular, the simulated DCC pions have significantly lower (pt) than the 
PYTHIA pions, so that the detection efficiencies ei and €2 are not the same. 
However, the values from PYTHIA are used here to illustrate the problems 
which will arise in DCC searches using these non-robust ratios.] 

The Ti^i obtained by counting numbers of charged tracks and photons 
aimed into the acceptance by PYTHIA are within two standard deviations 
of 1.00. This suggests that something close to a binomial charged- neutral 
distribution is used by PYTHIA. Detector effects are included by running 
the events through GEANT. The fact that the lower-order (statistically sig- 
nificant) Tj 1 did not differ by more than about 10% (and are in fact within 
2a) of the values from the PYTHIA input justifies the claim that the robust 
observables are indeed insensitive to detection efficiencies, and also that cor- 
relations of efficiencies with multiplicities and momentum do not greatly alter 
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the robustness of these variables. The ratios rj^i obtained from running the 
DCC events through GEANT are somewhat higher than the values predicted 
for the 1/(2a/7) distribution; however, they are clearly distinguishable from 
the values for generic production. 

The choice of parameters used in the DCC generator is somewhat opti- 
mistic; the DCC domain is aimed directly at the center of the acceptance, 
and the DCC pions have a rather large {pt)- As is shown in Fig. 4.17 , very 
low-pT photons are less likely to convert, and the efficiency for detecting those 
which do convert is only about 40% for px ^ 50 MeV. Charged pions with 
low Pt are also more difficult to find because they tend to stop or multiple 
scatter more in the lead; some scatter so much that the tail is not associated 
with the head by the vertexer, and may be found as a fake photon conver- 
sion. Therefore, we varied the {px) of the pions in the DCC generator. The 
resulting moments and rjj are given in Table The parameters in model 
A are those mentioned previously. Models B and C have lower {pt) , 50 MeV 
and 25 MeV, respectively. In these models, the ratio of mean energy density 
of the DCC pions to that of generic pions is not changed [if) = 1), which 
leads to larger numbers of pions in the domain since (iV^) oc ip/pr- In order 
to keep the mean number of pions constant, we also varied ip with pt in 
models D and E. For DCC domains with very low-pT pions, the values of the 
robust observables are closer to what we expect for a binomial distribution, 
so that distinguishing these domains from generically-produced pions using 
these observables becomes much more difficult. 

Of course, we do not expect to observe events consisting of only a DCC 
domain aimed into the acceptance. Possible scenarios for mixing DCC and 
generic multiparticle production are discussed in Ref. For example, 

any given event could be due to either DCC or generic production, but not 
both (exclusive production). Perhaps more realistically, the occurrence of 
DCC in an event could be independent of the generically-produced pions 
(independent production), or the amount of DCC production could depend 
on the amount of generic production (associated production). An example 
of the latter type is the Baked Alaska model |T^, which has the number of 
DCC pions given by Ndcc ~ Ngen^^^. 

In order to study the effect of an admixture of DCC with generic events 
where the amount of DCC produced is independent of the amount of generic 
production, we added DCC domains from the DCC-generator/GEANT to 
various fractions of (random) PYTHIA/GEANT events. The effect on the 
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Tj^i is shown in Table |6.3| . (Note that a shghtly older version of the vertexer 
was used for this study, which accounts for the discrepancy in ri^i here with 
no DCC added and in Table |0|.) The values for rj i fall off faster than 
linearly with fraction of DCC. For a fraction of 1, i.e. when a DCC domain 
is added to every event, the ratios are higher than for DCC alone, but are 
still easily distinguishable from the generic values. 



6.1.2 Characteristics of events from lead-in runs 

Ten of the lead-in runs (1096, 1099, 1103, 1109, 1110, 1125, 1126, 1127, 
1137, and 1139), totaling about 1.5 million events, were used in the following 
analysis. The running conditions were very clean (e.g. low luminosity) and 
therefore the diffractive tags had very little contamination from beam-gas 
interactions. The frequency distribution of events with given numbers of 
charged tracks and photons is determined (Table |6^) and used to calculate 
the factorial moments and Vij. Table |6.5| gives the values for some of these 
variables. The mean number of charged tracks found per event is about 0.5 
and of converted photons is about 0.2. The lower-order i are close to what 
is expected for a binomial distribution (r^ i = 1). The values for ri i and ri 2 
are within two standard deviations of the PYTHIA results. The higher-order 
ratios are weighted towards bins of M{nchi n^) which are statistically limited, 
and therefore the deviations from unity are not very significant. In any case, 
the ratios are not smaller than one as would be expected for a contribution 
from DCC. Therefore the events appear to be consistent with production by 
only generic mechanisms. 



6.1.3 Characteristics of events with diffractive tags 



From the 10 runs used in Sec. |6.1.2| , 21412 events have a ktag from the scin- 



tillator which detected showers from interactions of diffractive anti-protons 
in the kicker magnets. The mean numbers of charged particles and photons 
are lower for ktag events, as would be expected for diffractive events, where 
a large fraction of the total energy is carried away by the beam remnant, 
and the charged-charged and charged-gamma correlations are correspond- 
ingly lower. Table ^]5| gives the values of the r, ,,- for the ktag events. 

The upstream hadronic calorimeters at 2; ~ —25 m are used to tag events 
with diffractive anti-protons with xp ^ 0.5 and anti-neutrons. Differences 
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related to isospin exchange in diffractive events might be apparent in com- 
parisons between events with an li and those with a p. Figure 6^ shows 
histograms of the ADC readout from the li and p calorimeters. Events with 
ADC > 400 (in order to cut out background from products of showers in the 



magnets) were used to calculate the moments and rj given in Table 



The mean number of particles found is higher than that in events with the 
ktag, but lower than in the total sample of events, and is lower for the tag 
on leading li's than for p's with half the beam momentum, consistent with 
energy conservation. 

The for diffractive-tagged events do not differ by more than two stan- 
dard deviations from the values for the total sample. Therefore, we conclude 
that there is no evidence for more DCC production in events with diffractive 
tags, offering no support to the conjecture that Centauros are related to DCC 
and are diffractive in nature. 



6.1.4 Characteristics of events with a pbar multiphcity 
tag 

Since DCC may be more often present in events with large multiplicity, we 
would like to find a measure of the multiplicity independent of that in the 
small acceptance. The multiplicity in the scintillator on the downstream 
anti-proton side of the collision at z = —81 in ("pbar counters") may be 
correlated with the total multiplicity, and therefore could be such a measure. 

In GEANT, where the number of charged tracks hitting the pbar counters 
is known, a plot of the total energy deposited in the pbar counters against 



the pbar multiplicity (Fig. |6.3|) shows the correlation between energy and 
multiplicity which we expect from minimum-ionizing particles. The pbar 
multiplicity in the real data can thus be found from the mip peaks in the 
ADC readout of the counters (Fig. |6^ ). 

In this way, the pbar multiplicities of events in runs 1099, 1103, 1109, 
1125, 1126, and 1127 were determined, and events were grouped in bins of 
pbar multiplicity, such that each bin contained 10% of the events. The multi- 
plicity in the pbar counters is indeed correlated with that in the acceptance. 
Table |6.6| shows an increase in (rich) and (n^) with increasing pbar multi- 
plicity. However, the rj_i do not appear to vary with multiplicity, and show 
no sign of an increased presence of DCC for higher multiplicity events. The 
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rial's are consistent with each other and with that for all events. 



6.2 Multiparticle analysis 

Plots of dNch/dt] and dN^/di], uncorrected for detection and trigger efficien- 
cies, are shown in Fig. |6.5| and Fig. respectively. The solid line is the 
result from approximately 1.5 x 10^ events from the 10 runs used in previous 
analysis. The dotted line is the result from the PYTHIA/GEANT simu- 
lation, and the dashed line is from the PYTHIA events which passed the 
GEANT trigger (about 1.5 x 10^ events). The charged multiplicity found 
by the vertexer for GEANT events is higher than that from the PYTHIA 
input by about 0.15, or 3%. The number of photons found for GEANT is 
about 42% of the PYTHIA input. The charged multiplicity from the data is 
lower than that from GEANT and from PYTHIA (about 90% of the GEANT 
value). The photon multiplicity is approximately 92% of that in GEANT. 

Since multiplicity is correlated with NHITS, and the contribution to the 
NHITS from background is greater in the real data than in GEANT, it is not 
unlikely that the cut of NHITS < 600 would cut more high-multiplicity events 
from the data sample than from the GEANT. This conjecture is supported by 
the results in Sec. [4.3.5| , which showed a correlation between pipe shower and 



multiplicity. However, we do not necessarily expect the data to agree with 
the simulations, since the PYTHIA input is not based on real measurements 
at the cm energy and 1] range of the MiniMax experiment. 



6.3 Low pt photons 

There has been controversial evidence suggesting that the number of low-p^ 
photons produced is larger than what is expected from hadronic decays and 
QED inner bremsstrahlung, in particular, that there is an excess of photons 
with Pt ^ 20 MeV [B^ . We have attempted to study this by observing pho- 
ton conversions in the 1/4 in (0.07 Xq) -thick Al window at z = 120 in, and 
measuring the momentum of the conversion products either by their bending 
in the scintillator (0.03 Xg-thick at z = 157 in) due to multiple scattering or 
by tracking them into the calorimeter and determining the energy deposited 
in the hit cells. The lead-out runs were used for this so that the conversion 
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tracks would not stop or scatter and lose energy in the lead. Conversions 
in the window, rather than in the scintillator, for example, were chosen be- 
cause the increased resolution in v in the front chambers was desirable for 
determining the momentum of the converted photon. 

First, it was necessary to determine which tracks were due to window 
conversions. To do this, we plotted the z of closest approach between pairs 
of tracks. The GEANT lead-out data shows a wide peak at the collision 
point, and a smaller peak at the window, z ^ 120 in, as can be seen in Fig. 
|6.71 . For the real data, the peak at the collision point and at the window 
are fairly clear, but there is also a smaller peak at intermediate z. This is 
apparently due to fake tracks made from u wires borrowed from several real 
tracks, such that the resulting fakes tended to vertex with the real tracks at a 
location in z between the collision point and the first chamber (which is just 
behind the window) . These fake vertices are removed by requiring the heads 
to point such that they will hit the lead, and that they hit all eight of the 
front chambers. The distribution of the z of closest approach for remaining 



tracks is shown in Fig. 5.8 



A window conversion is defined as a pair of tracks with a z at closest 
approach between 115 and 150in, and with a reduced of the match between 
the tracks at the window (done in the same way as the usual vertexing at the 
lead) of xV2 < 2.5 (see Fig. ^ for the GEANT distribution and Fig. |]8 



for that of the data). Also, in order to be able to determine the momentum 
of the conversion tracks by bending in the scintillator or by following them 
into the calorimeter, each head is required to have exactly one tail which 
matched to the head with < 5. 

Both conversion tracks are found by the tracker in about 26% of the 
window conversions in GEANT which send both tracks into the acceptance. 
Only one track is recorded 70% of the time, usually because the tracks are so 
close together that they are grouped by the tracker into a single track. This 
is apparent in Fig. |6.9| , which shows the opening angle of window conversion 
tracks in GEANT and of the pairs of tracks found as window conversions. 
Conversions with smaller opening angles are much less likely to be found. 

Of the pairs of tracks found as window conversions, approximately 90% 
have one tail matched to each head. As is also shown in Fig. |6]^, the efficiency 
for finding vertices from the low-pT photons that convert in the window is 
fairly high (because such conversions tend to have larger opening angles). 
Unfortunately, the low-p^ conversion tracks scatter more in the scintillator 
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than high energy particles, so that the efficiency for matching the tails to 
the heads is relatively low, preventing a determination of the momentum of 
these tracks. 

The momentum of the conversion tracks is determined from the angles 
between the heads and tails at the scintillator using the relationship between 
the momentum of a particle and the bending angle of the track due to mul- 
tiple scattering in a material pBl 



^ 13.6 MeV [IF 

Pep 'ixo' ^^-^^ 

where x/Xq is the thickness of the material in radiation lengths, p is the 
momentum and (3c the velocity of the particle, z is its charged number, and 
^0 is the bending angle in the plane of the incident and bent track. The 
momentum of the photons which converted in the window is taken to be the 
sum of the momenta of the two conversion tracks. The pt is determined from 
the total momentum and the angle of the trajectory of the photon, which 
must have traveled from the collision point to the vertex of the conversion 
tracks. For GEANT conversions, where the photon momentum is known, the 
Pt as determined by bending in the scintillator tends to be lower than the 
actual Pt (Fig. |6.1CI| ). Other attempts to find a measure of the momentum, 
using, for example, the opening angle, were unsuccessful. 

Following the conversion tracks into the calorimeter, their energy is de- 
termined to be simply the energy deposited in the cell(s) into which the 
tracks are aimed, without using a clustering algorithm. The pt found with 
this method is surprisingly accurate (Fig. |6.1U| ). Unfortunately, the actual 



calorimeter was never successfully calibrated, and therefore this method can 
not be reliably used on the data. 

The Pt distributions of "photons" from pairs of tracks which are found 
as window conversions in the data are shown in Fig. |6.11|, although these 



distributions are not believed to be very accurate. Instead of comparing pt 
distributions from GEANT and from the data, a comparison of the number 
of window conversions integrated over all pt was used to look for evidence 
for an excess of low-pr photons. 

First, the origins of pairs of tracks found as GEANT window conversions 
are determined. The relative contribution from sources other than actual 
window conversions should be low, and hopefully from sources of similar 
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magnitude in the GEANT and in the data. In a sample of about 10^ GEANT 
events, 225 pairs of tracks are identified as window conversions. About half 
of these really are window conversions (106 conversions in the window and 4 
conversions in the first chamber at z = 123 in). The next largest contribution 
is 19% from pairs of charged tracks from the collision point which are close 
together at the window, including 13 dalitz pairs. The large uncertainty in 
the point of closest approach due to the nearly identical slope of such tracks 
allows the possibility of their appearing to vertex at the window. Strong 
interactions in the window, neutral particle decays (e.g. K°, A°), pairs of 
a primary charged track from the collision and a secondary charged track, 
and pairs of a charged track and a fake each contribute between 5 and 7% 
of the pairs found as window conversions. Five events have conversions in 
the window of secondary photons. The remaining few events involve pairs of 
secondary tracks from interactions in other material. 

Almost 10^ events from lead-out runs (1089, 1093, 1104, 1108, 1111, 1123, 
1124, 1129, 1132) were used to search for window conversions in the data. 
The ratio of the number of window conversions found in the data to that in 
GEANT is 2121/339 = 6.3. The ratio of events in the data sample to that 
in the GEANT sample is 923238/155737 = 5.9, so that there appears to be 
an excess of window conversion in the data of about 6%. The uncertainties 
involved in this analysis do not allow a precise statement either supporting 
or contradicting the claim of Ref . [^] . 
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(r(i)-r(mean))/sigr(i) 



Figure 6.1: Histogram of (r(i) — {r))/ar{i) for groups {i) of 25000 events and 
Gaussian fit. 
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Figure 6.2: ADC values of the hadronic calorimeters which see anti-neutrons 
[ADC(40)] and anti-protons [ADC(42)]. 



130 





obar.hbook 








Dbar.hbook 




600C 






ID 101 
Entries 5071 7 
Meon 6.132 
RMS 3.799 


1 000 








ID 103 
Entries 50717 
Meon 43.68 
RMS 26.10 




5000 


















4000 








800 












3000 








600 












2000 








400 












1 000 








Z u w 















t 




1 1 1 1 1 1 1 1 









1 1 1 1 1 1 1 


1 1 1 1 1 1 







10 20 


30 





20 40 60 80 100 


50 


Ncharged at pbar 
Dbar.hbook 


50 




Epbar 

Dbar.hbook 




















45 








45 












40 








40 












35 








35 












30 








30 












25 








25 












20 








20 












15 








15 












10 








10 












5 








5 




























■y, , 1 , , , , 1 , 






2.5 5 


7.5 10 







2.5 5 


7.5 10 






Epbor vs Ncharged ot pbor 








Epbor vs Nchorged at pbar 



Figure 6.3: Energy (in MeV) deposited in the pbar counters vs number of 
charged particles hitting these counters 
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Figure 6.4: PED_SUB_ADC values (in units of mip energies or, equivalently, 
number of charged particles) of the pbar counters 
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Figure 6.5: Pseudorapidity distribution dN/dr] vs rj for charged tracks in 
the data (solid), GEANT (dotted), and PYTHIA (dashed), uncorrected for 
efficiencies. 
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Figure 6.6: Pseudorapidity distribution dN/dr] vs rj for photons in the data 
(sohd), GEANT (dotted), and PYTHIA (dashed), uncorrected for efficien- 
cies. 
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Figure 6.7: Histogram of the z of closest approach between all heads, and of 
the of the matching of those pairs in GEANT. 
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Figure 6.8: Histogram of the z of closest approach between all heads, and of 
the of the matching of those pairs in the data. 
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Figure 6.9: Opening angle between pairs found as window conversion tracks 
in GEANT (solid) and between all GEANT conversion tracks (dashed), px of 
all primary photons in GEANT, and efficiencies for finding window conver- 
sions as a function of pt- In the lower-left plot, the solid line represents the 
fraction of primary photons which converted in the window, the dashed line 
is the fraction which also produced two conversion tracks in the acceptance. 
In the lower right plot, the solid line is the fraction of window conversions 
with both tracks in the acceptance which were found as a pair of heads, and 
the dashed line is the fraction which also had one tail per head. 
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Figure 6.10: Contour plots of the pt found by bending in the scintillator 

and by using the calorimeter vs the actual pt of the converted photon in 
GEANT, and corresponding px distributions (solid is calculated pr, dashed 
is actual). 
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Figure 6.11: Distributions of px found by bending in the scintillator and by 
using the calorimeter for converted photons in the data. 
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Table 6.1: Robust observables r^j for generic events simulated by PYTHIA 
and pure-DCC events simulated by the DCC generator, along with predic- 
tions for binomial and l/{2^/f) distributions. 
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Table 6.2: Results for varying parameters in the DCC generator. 
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Table 6.3: The effect on the i of an admixture of DCC and generic 
(PYTHIA) events. DCC domains from the DCC generator/GEANT are 
added to various fractions of random PYTHIA/GEANT events. The first 
column represents the fraction of events in which a DCC is overlaying a 
generic event. A DCC fraction of 1 means that DCC has been added to 
every event, not that the events are pure DCC as in Table 6.1. 
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Table 6.4: Number of events with a given rich, ti. 
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1 052 ± 036 

X . W ^ —1— KJ • kJKJkJ 


1 051 ± 031 


1 076 ± 024 


r(2 I) 


1 0320 + 00Q7 


1 1 34 + n 1 07 


1 042 + 083 




r(3 1) 


1 0563 ± 0251 


1 302 ± 281 


1 220 ± 211 

X • ^ ^ W 1 V7 • ^ X X 


1 180 ± 173 

X . X w w — 1 — w . X i o 


r(4, 1) 


1.1183 ±0.0612 


1.802 ± 0.637 


2.102 ±0.599 


1.370 ±0.413 


r(5,l) 


1.3382 ±0.1436 






1.721 ±0.911 


r(6,l) 


2.0758 ±0.4159 








r(7, 1) 


5.0105 ± 5.6885 








r(0.2) 


1.5793 ± 0.0092 


1.647 ±0.090 


1.683 ±0.077 


1.658 ±0.059 


r(l,2) 


1.5632 ±0.0204 


1.623 ±0.217 


1.700 ±0.176 


1.675 ±0.152 


r(2,2) 


1.5538 ±0.0497 


1.914 ±0.616 


1.908 ± 0.474 


1.933 ±0.441 


r(0,3) 


3.3711 ±0.0759 


3.378 ±0.627 


3.359 ± 0.545 


3.333 ±0.390 


r(l,3) 


3.1115 ±0.1210 


3.335 ± 1.165 


3.196 ±0.895 


3.367 ±0.958 


r(0,4) 


8.9729 ± 0.6201 


6.675 ± 2.934 


7.266 ± 2.980 


6.236 ± 1.780 



Table 6.5: Values of for all lead-in events and for those with diffractive 
tags. 
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bin 


{rich) 




'^1,1 


r2,i 




1 


0.364 ± 0.002 


0.145 ± 0.002 


1.02 ±0.02 





97 ±0.06 


0.87 ±0.14 


2 


0.401 ± 0.002 


0.165 ± 0.002 


1.00 ±0.02 


1 


00 ±0.05 


1.04 ±0.14 


3 


0.426 ± 0.002 


0.173 ±0.002 


1.03 ±0.02 


1 


11 ±0.05 


1.28 ±0.13 


4 


0.449 ± 0.002 


0.183 ±0.002 


1.00 ±0.02 


1 


05 ±0.05 


1.10 ±0.13 


5 


0.468 ± 0.003 


0.190 ±0.002 


1.03 ±0.02 


1 


10 ±0.05 


1.17 ±0.14 


6 


0.493 ±0.003 


0.195 ±0.002 


1.02 ±0.02 


1 


01 ± 0.04 


1.05 ±0.11 


7 


0.503 ± 0.003 


0.202 ± 0.002 


1.04 ±0.02 


1 


03 ± 0.05 


0.98 ±0.12 


8 


0.531 ±0.003 


0.215 ±0.002 


1.01 ±0.01 


1 


00 ± 0.04 


0.97 ±0.09 


9 


0.555 ±0.003 


0.221 ±0.002 


1.03 ±0.01 


1 


06 ± 0.04 


1.09 ±0.11 


10 


0.603 ± 0.003 


0.237 ±0.002 


1.03 ±0.01 


1 


02 ± 0.03 


1.01 ±0.07 



Table 6.6: Mean multiplicities and robust observables for bins of pbar mul- 
tiplicities each containing 10% of the events. The pbar multiplicity increases 
with increasing bin number 1-10. 



Chapter 7 
Conclusions 



MiniMax ran successfully with many different detector configurations for a 
period of about two years. During this time, much was learned which should 
prove useful for the operation and design of future detectors in the forward 
region. More than 1.5 x 10^ events with lead in and almost 10^ events with 
lead out from clean runs when the detector was running properly have been 
analyzed for the work presented here. 

The GEANT simulation does not accurately represent the amount of 
background seen in the wire chambers from interactions in material surround- 
ing the detector. Even with all of the material included in the simulation, 
the mean number of hits in the MWPC's is low by about a factor of two. 
However, the PYTHIA/GEANT output is useful for setting upper limits on 
how well the data analysis tools are able to reconstruct events. Based on 
the efficiencies and numbers of fakes discussed in Sec. [4.3| , the tracker and 
vertexer appear to work quite well. 

The simulations also show that the robust observables of Chapter 5 are 
insensitive to detector effects and physics complications, as advertised. The 
PYTHIA/GEANT simulation, which includes a photon conversion efficiency 
of slightly better than 50%, an efficiency of 80% for detecting converted 
photons, resonance production, simulations of detector effects, and many 
other features, matches the predictions of a simple binomial model. At least 
for DCC modeled as in Sec. |3.3| , the robust observables should be useful in 
distinguishing DCC from generic production. 

We have found no evidence for DCC in the total sample of events or in 
diffractive events. Limits on the amount of DCC which could be present 
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without being detected depend on the models of DCC (e.g. the Pt of the 
DCC pions) and of combining DCC with generic production. 

Raw measurements of dNch/dr) and dN^/drj have been made in a previ- 
ously unexplored region of phase space, from which actual distributions can 
be derived, taking into account detection efficiencies and fakes, and trigger 
efficiencies. Further work on modeling fakes in the data is necessary before 
this can be done. 



Appendix A 

Uncertainty Calculations 



The normalized factorial moments are defined as 

{ndnc - 1) . . . (ric - i + l)ng{ng - 1) . . . {ug - I + 1)) 



{ucY {rig) 



(A.l) 



where (O) = ^ Enc,% C'A/'K, n^), N = Enc,ng^{nc,ng), and A/'K, n^) is 
the number of events with ric charged tracks and Ug gammas. 

We assume that = N{nc,ng). The uncertainty in Fj^i is then 

given by 



{[uc . . . (nc - J + l)ng ...{ug - I + 1)]^) ^yp, 2 



{ric 



-2j 



. {ndric ...{ric- J + 1)% ■ ■ ■ {rig - I + 1)]) 



o, KK ■ ■ ■ K - J + l)ng ...{ng-l + 1)]) 

— Z( m -T i J 



(A.2) 
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In order to find the correlated uncertainty for any two F's we define 

X = Fj,i + Fi,m (A.3) 
which has uncertainty given by 

(^l = 4,, + 4,^ + 2ai.,,F,^ (A.4) 
Writing the F's in terms of the M{nc, Ug), we get 
"<...«-j + lK...(n' -/+1) 



+ + 



N 



^N{nc) N{ng) 



^ + l)^g ...{n'g-m + l) 

N{nS{ngr 



-I——, — : — m-^ 



n 2 



iV(ne) N{ng)_ 
and we can pick out the correlation terms to find 



(A.5) 



1 

TV 



{[uc ...{nc-j + l)ng . ..{ng-l + l)][nc. . .{nc-i + 1)% ...{ng-m + 1)]) 



. (njric . . . (ric - j + 1)^^ . . . (n^ - / + 1)]) 



{ugY 



—m 



{ng[nc ...{nc-j + l)ng . . . {ug - I + 1)]) 



{ucY {rig) 



. (ncK ■ ■ ■ (wc - i + l)ng . . .{ug -m + 1)]) 



{UcY {ng) 
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{ng[nc ...{nc-i + l)ng ...{ng-m + 1)]) 



{Uc) 



ill 11 I 

+{jm + tl) ^ " - (j + / - l)(i + m - l)F,yF,, 



(A.6) 



The generalized r defined by 



r,i = ^ (A.7) 



has uncertainty given by 



Appendix B 
Track fitting 



B.l Track fitter: 

SUBROUTINE fit_init 
IMPLICIT NONE 

INTEGER I,NN 

REAL sigsq,wiresp ,pi ,uajig, sum, cu, su 
REAL U0(24) ,V0(24) ,aiig(24) ,Z(24) 

common/ cham_param/NN , sigsq , wiresp , pi , uang , sum , cu , su , UO , VO , ang , Z 
C position of the center of wire an.d angle of rotation in (x,y,z) 



real zz (24) /123 . 69 , 126 . 57 , 129 . 69 , 132 . 57 , 
+ 135.69,138.69,141.69,144.69, 
+ 167.09,169.89,172.84,176.09, 
+ 178.84,181.84,185.09,189.65, 
+ 194.29,197.29,200.29,203.29, 
+ 206.16,209.16,212.41,215.41/ 

real xxO (24) /2 . 806 , 10 . 1 17 , . 853 , 11 . 792 , 
+ 9.997,0.818,10.076,1.252, 
+ 1.692,10.090,10.970,1.535, 
+ 11.199,12.360,11.424,3.366, 
+ 11.602,11.327,2.814,12.648, 
+ 12.056,3.560,3.140,2.681/ 

real yyO (24) /I 1 . 653 , 1 1 . 032 , 10 . 870 , 9 . 255 , 
+ 10.910,9.193,11.700,1.813, 
+ 2.568,12.305,11.702,3.508, 



code 

fit.f 
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+ 11. 781, 10. 771, 11-812, 1-695, 

+ 12.036,12.524,2.963,11.415, 

+ 12.114,1.613,3.048,3.748/ 

real aangle (24) /-I . 0315 , -2 . 3422 , -0 . 7575 , -2 . 5953 , 

+ -2.3318,-0.5044,-2.0857,0.8029, 

+ 0.7959,-2.1084,-2.3370,0.5794, 

+ -2 . 3370 , -2 . 5063 , -2 . 3405 , . 9669 , 

+ -2 . 3370 ,-2.2201,0. 8029 , -2 . 4574 , 

+ -2.3527,1.0420,0.7837,0.6266/ 



sigsq=l./12. luncer in location of track (wire spacings) for each hit 

wiresp=0.1 !wire spacing in inches 

NN=4 

pi=3. 141592654 

uang=45.43*pi/180. ! taken as meein eingle of "u chambers" 

cu=cos(uang) 

su=sin(uang) 

do i=l,24 

UO(i)=xxO(i)*cu+yyO(i)*su !u,v of center of wire 

VO(i)=-xxO(i)*su+yyO(i)*cu 

Z(i)=zz(i) 

if (aangle (i) . It . 2*pi) aangle(i)=aangle(i)+2*pi 

axLg ( i ) =aaiigle ( i ) -uang 

ang(i)=ang(i)-pi/2. 

enddo 

RETURN 
END 



SUBROUTINE fit_wires 
IMPLICIT NONE 
INTEGER I,J,K,L,nijm 

INTEGER icham,ichlast,ncham,iwire,nwires(24) ,iichaiii(24) 
REAL wire 

REAL awires(24,128) 

REAL*8 X(4,4) ,XX(100,4) ,Y(4) ,YY(100) ,a(4) ,XINV(4,4) 
REAL*8 INDEX (4), D 
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INTEGER NN 

REAL s igsq , wiresp , pi , uEing , sum , cu , su 

REAL U0(24) ,V0(24) ,ang(24) ,Z(24) 

common/ cham_pareim/NN , sigsq , wiresp , pi , uang , sum , cu , su , UO , VO , ang , Z 

INTEGER HLIM.HBSIZE 
PARAMETER (HBSIZE=500000) 
COMMON/PAWC/HLIM (HBSIZE) 

include ' wires. inc' 
include 'fit. inc' 

include ' track_dat . inc ' 
include 'dst.inc' 

integer ib(24) ,dmin,newwire 

include 'EVENT. INC !to find missing hit in front-8 tracks 

DO I=l,ntrack 
100 continue 

C will fit track parameters 'a' in matrix eq Xa=Y 

DO J=l,4 

DO K=l,4 

X(J,K)=0 

ENDDO 

Y(J)=0 
ENDDO 

C will need to know which chambers and which wires in each chamber 
C were hit to get chisq 

nchajii=0 

ichlast=0 

if (nnwire(i) .gt . 100) write(6,*) 'nnwire=' ,nnwire(i) 
DO J=l ,nnwire(i) !go through hit wires from dst 

icham=int ((wiresCl, J)-l)/128. )+l 
iwire=wires (I , J) - ( icham-1 ) * 128 
IF (icham.ne.ichlast) THEN 
ncham=ncham+ 1 
iicham(ncham) =icham 
nwires (ncham) =1 
ELSE 

nwires (ncham) =nwires (ncham) +1 
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ENDIF 

awires (ncham , nwires (ncham) ) =iwire 
ichlast=icham 

C fit XX(i, j)a(j)=YY(i) for all chambers i to get a(j) 
C min chisq -> Y=Xa, a=(Xinv)Y 

C (take Scime uncertainty for all hits siqsq=l/12 (wiresp~2)) 

XX(J, l)=-sin(ang(icham)) 
XX( J, 2)=-Z(ichciin) *sin(ang(icham) ) 
XX(J,3)=cos(ang(ichain) ) 
XX(J,4)=Z(ichain)*cos(ang(ichajn) ) 

YY( J)=-UO(ichciin) *sin(ang(ichciin) )+VO(ichain)*cos(ang(ichain) ) + 
+ iwire*wiresp 
DO L=l,4 
DO K=l,4 

X(L,K)=X(L,K)+XX(J,L)*XX(J,K) 
ENDDD 

Y (L) =Y (L) +XX ( J , L) *YY ( J) 
EMDDO 
999 continue 
ENDDO 

C invert X 

DO J=l,4 
DO K=l,4 

XINV(J,K)=0 
ENDDO 

XIMVCJ, J)=l 

ENDDD 

C Numerical Recipes subroutines to invert a matrix 
CALL LUDCMP(X,NN, INDEX, D) 
DO K=l,4 

CALL LUBKSB(X,NN, INDEX, XINVd.K)) 
ENDDO 

DO J=l,4 
a(J)=0 
DO K=l,4 

a(J)=a(J)+Y(K)*XINV(J,K) !get track parameters 
ENDDO 
ENDDO 



C fill fit.inc 
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a_u(I)=a(l) 
b_u(I)=a(2) 
a_v(I)=a(3) 
b_v(I)=a(4) 

C get chisq (if more than one wire in a chamber is part of the track, 

C use residuals from each wire 

C total number of wires hit = num) 

chisq(I)=0 

num=0 

DO J=l,nchain 
K=iicheim( J) 

wire=(-(a_u(I)+b_u(I)*Z(K)-UO(K))*sin(ang(K)) 
+ + (a_v (I) +b_v (I ) *Z (K) -VO (K) ) *cos (ang (K) ) ) /wiresp 
DO L=l,nwires(J) 
num=num+ 1 

chisq(I) =chisq(I)+ (wire-awires ( J ,L) ) **2 
ENDDO 
ENDDO 

chisq(I)=chisq(I)/(sigsq*(num-4)) 

C get covariance matrix 
DO J=l,4 
DO K=l,4 

CO var ( I , J , K) = ( s igsq*0 . 01) *XINV ( J , K) 
ENDDO 
ENDDO 

enddo 
endif 

ENDDO 

RETURN 
END 

C 

SUBROUTINE LUDCMP(A,N1,INDX,D) !LU decomposition 

IMPLICIT NONE 

integer n,i, j ,k,imax,Nl 

real*8 nmax,tiny,aamax,sum,dum Ineed double precision to get XinvX=l 



155 



PARAMETER (NMAX=46 , TINY=1 . OE-20) 
REAL*8 NP,D,A(4,4) ,IMDX(4) ,VV(4) 

N=N1 
D=l. 

DO 12 1=1, M 
AAMAX=0 . 
DO 11 J=1,M 

IF (ABS(A(I, J)) .GT.AAMAX) AAMAX=ABS(A(I, J)) 

11 CONTINUE 

IF (AAMAX.EQ.O.) PAUSE 'Singular matrix. ' 
VV(I)=1./AAMAX 

12 CONTINUE 

DO 19 J=1,N 
IF (J.GT.l) THEN 
DO 14 1=1, J-1 
SUM=A(I,J) 
IF (I.GT.l)THEN 
DO 13 K=1,I-1 

SUM=SUM- A (I,K)*A(K,J) 

13 CONTINUE 
A(I,J)=SUM 

ENDIF 

14 CONTINUE 
ENDIF 
AAMAX=0 . 

DO 16 I=J,N 
SUM=A(I,J) 
IF (J.GT.l)THEN 
DO 15 K=1,J-1 

SUM=SUM-A ( I , K) *A (K , J) 

15 CONTINUE 
A(I,J)=SUM 

ENDIF 

DUM=VV(I)*ABS(SUM) 

IF (DUM.GE.AAMAX) THEN 

IMAX=I 

AAMAX=DUM 
ENDIF 

16 CONTINUE 

IF (J.NE.IMAX)THEN 
DO 17 K=1,N 
DUM=A(IMAX,K) 
A(IMAX,K)=A(J,K) 
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A(J,K)=DUM 

17 CONTINUE 

D=-D 

VV(IMAX)=VV(J) 
ENDIF 

INDX(J)=IMAX 

IF(A(J,J) .EC!.0.)A(J,J)=TINY 
IF(J.NE.N)THEN 

DUM=1./A(J, J) 

DO 18 I=J+1,N 

A(I,J)=A(I,J)*DUM 

18 CONTINUE 
ENDIF 

19 CONTINUE 
RETURN 
END 



SUBROUTINE LUBKSB(A,N1,INDX,B) Iback substitution 

IMPLICIT NONE 
integer n,ii,i,ll, j ,N1 
REAL*8 sum,NP,A(4,4) ,INDX(4) ,B(4) 
11=0 
N=N1 

DO 12 1=1, N 
LL=INDX(I) 
SUM=B(LL) 
B(LL)=B(I) 
IF (II.NE.O)THEN 
DO 11 J=II,I-1 

SUM=SUM-A(I,J)*B(J) 

11 CONTINUE 

ELSE IF (SUM.NE.O.) THEN 

II=I 
ENDIF 
B(I)=SUM 

12 CONTINUE 

DO 14 I=N,1,-1 
SUM=B(I) 
IF(I.LT.N)THEN 
DO 13 J=I+1,N 

SUM=SUM-A(I,J)*B(J) 

13 CONTINUE 
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ENDIF 

B(I)=SUM/A(I,I) 

14 CONTINUE 
RETURN 
END 



C 

C wires. inc (filled by eg vertexer.f) 

integer ntrack , nnwire (200) , wires (200 ,100) 

common/wiresf orf it/ntrack , nnwire , wires 

C (must be filled by user unless using the vertexer) 

C ntrack=# tracks 

C nnwire (i)=# wires in track i 

C wires (i,j)=jth wire of ith track in usual format (icham-1) *128+iwire 



C 

C fit. inc (filled by fit.f) 

REAL a_u(200) ,b_u(200) ,a_v(200) ,b_v(200) , chisq(200) 

REAL C0VAR(200,4,4) 

integer current_evnum 

common/track_parain/a_u , b_u , a_v , b_v , chi sq , covar , current _evnum 

C (determined by fit.f) 

C u=a_u+b_u*z, v=a_v+b_v*z 

C 

C for covaricince matrix, indices go like 
C 1 = a_u, 2 = b_u, 3 = a_v, 4 = b_v 
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B.2 Code for correcting tracks with double 
hits in the non-u chambers: refit. f 

SUBROUTINE refit 
IMPLICIT NONE 
INTEGER I,J,K,L,M 

INTEGER icham, ichlast .ncham, iwire .nwires (24) , iich.ciin(24) 
INTEGER NN 

REAL sigsq , wiresp , pi , UEing , sum , cu , su 
REAL U0(24) ,V0(24) ,aiig(24) ,Z(24) 

common/ cham_param/NN , sigsq , wiresp , pi , uang , sum , cu , su , UO , VO , ang , Z 
integer itrack 

real av,bv,v(24) ,u(24) ,vlead,ucham(24) 

real sz,szz,sv,szv,cth(24) ,sth(24) ,denom,predwire, wire (18000, 24) 
real sigth.thgam.thtail 

integer nwirepl(24) , iwirepl(24, 10) ,nwire,icomb 

integer iw(8) ,il, 12, 13, 14, 15, 16,17,18 

Integer ncomb , ichi , ipoint 

real chi(18000) .point (18000) , chimin, pointmin 

integer ntr,hit2 
integer nchl,nch2 
common/newnch/nchl ,nch2 

include ' wires. inc' 
include 'fit. inc' 

include ' track_dat . inc ' 
Include 'dst.lnc' 

C get u chambers 
do j=l,24 
ucham(j)=0 
enddo 

ucham(2)=l 
ucham(5)=l 
ucham(8)=l 
ucham(9)=l 
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ucham(ll)=l 
uchain(13)=l 
uch.cim(15)=l 
ucham(17)=l 
uchain(19)=l 
ucliam(23)=l 
uchain(24)=l 

C will need orientation of non-u chambers to fit in v 
do k=l,24 

if (uchain(k) .eq.O) then 
sth(k)=sin(ang(k) ) 
cth (k) =cos (ang(k) ) 
endif 
enddo 

C check number of wires hit per plane 

C if 2 wires hit are separated by more thein 3 wires in amy pleine, refit 
do i=l,ntrack 
do k=l,24 

nwirepl(k)=0 
enddo 

do j=l ,nnwire(i) 

icham=int ( (wires (i , j ) -1 . ) /128 . )+l 

nwirepl (icham) =nwirepl (ichcun) +1 

iwirepl (icham , nwirepl (icham) ) =wires ( i , j ) 
enddo 
hit2=0 
do k=l,24 

if ( (nwirepl (k) . eq. 2) . ctnd . 
+ (abs(iwirepl (k, l)-iwirepl(k, 2) ) . gt . 3) ) then 

hit2=l 

endif 
enddo 

if (hit2.eq.O) goto 999 

C how many chambers have "double hits" and which wires are involved 
ncham=0 
do k=l,24 

nwirepl (k)=0 
enddo 
ichlast=0 
DO J=l ,nnwire(i) 

icham=int((wires(I, J)-l)/128.)+l 
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if (ucham(icham) . eq. 0) then 

if (icham.ne . ichlast) nchain=nchain+l 

iicham(ncham) =ichaim 

nwirepl (ncham) =nwirepl (ncham) +1 

iwirepl (ncham , nwirepl (ncham) ) =wires ( I , J) - ( icham-1 ) * 128 
ichlast=icham 
endif 
ENDDO 

C cycle thru all wires ! ! ! 

C and find all possible combinations with 1 wire per plane 

ncomb=0 

if (ncham. eq. 3) then 
do i 1=1 .nwirepl (1) 

iw(l)=il 
do 12=1, nwirepl (2) 

iw(2)=i2 
do 13=1 , nwirepl (3) 

iw(3)=i3 

ncomb=ncomb+ 1 

do j=l, ncham 

if (ncomb.le. 18000) wire(ncomb, j)=iwirepl(j ,iw(j)) 

enddo 
enddo 
enddo 
enddo 
endif 

if (ncham. eq. 4) then 
do 11=1 .nwirepl (1) 

iw(l)=il 
do 12=1. nwirepl (2) 

iw(2)=i2 
do 13=1 ,nwirepl(3) 

iw(3)=i3 
do 14=1. nwirepl (4) 

iw(4)=i4 

ncomb=ncomb+l 

do j=l. ncham 

if (ncomb.le. 18000) wire(ncomb, j)=iwirepl(j ,iw(j)) 

enddo 
enddo 
enddo 
enddo 
enddo 
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endif 

if (ncham. eq. 5) then 
do il=l ,nwirepl(l) 

iw(l)=il 
do i2=l ,nwirepl(2) 

iw(2)=i2 
do ±3=1 ,nwirepl(3) 

iw(3)=i3 
do i4=l ,nwirepl(4) 

iw(4)=i4 
do i5=l ,nwirepl(5) 

iw(5)=i5 

ncomb=ncomb+l 

do j=l, ncham 

if (ncomb.le. 18000) wire(ncomb, j)=iwirepl(j ,iw(j)) 

enddo 
enddo 
enddo 
enddo 
enddo 
enddo 
endif 

if (ncham. eq. 6) then 
do il=l ,nwirepl(l) 

iw(l)=il 
do i2=l ,nwirepl(2) 

iw(2)=i2 
do i3=l ,nwirepl(3) 

iw(3)=i3 
do i4=l ,nwirepl(4) 

iw(4)=i4 
do i5=l ,nwirepl(5) 

iw(5)=i5 
do i6=l ,nwirepl(6) 

iw(6)=i6 

ncomb=ncomb+l 

do j=l, ncham 

if (ncomb.le. 18000) wire(ncomb, j)=iwirepl(j ,iw(j)) 

enddo 
enddo 
enddo 
enddo 
enddo 
enddo 
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enddo 
endif 

if (ncham.eq.7) then 
do il=l ,nwirepl(l) 

iw(l)=il 
do i2=l ,nwirepl(2) 

iw(2)=i2 
do i3=l ,nwirepl(3) 

iw(3)=i3 
do i4=l ,nwirepl(4) 

iw(4)=i4 
do i5=l ,nwirepl(5) 

iw(5)=i5 
do i6=l ,nwirepl(6) 

iw(6)=i6 
do i7=l ,nwirepl(7) 

iw(7)=i7 

ncomb=ncomb+l 

do j=l,ncham 

if (ncomb.le. 18000) wire(ncomb, j)=iwirepl(j ,iw(j)) 

enddo 
enddo 
enddo 
enddo 
enddo 
enddo 
enddo 
enddo 
endif 

if (nchcim.eq.8) then 
do il=l ,nwirepl(l) 

iw(l)=il 
do i2=l ,nwirepl(2) 

iw(2)=i2 
do i3=l ,nwirepl(3) 

iw(3)=i3 
do i4=l ,nwirepl(4) 

iw(4)=i4 
do i5=l ,nwirepl(5) 

iw(5)=i5 
do i6=l ,nwirepl(6) 

iw(6)=i6 
do i7=l ,nwirepl(7) 

iw(7)=i7 
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do 18=1 ,nwirepl(8) 
iw(8)=i8 
ncomb=ncomb+l 
do j=l,ncham 

if (ncomb.le. 18000) wire(ncomb, j)=iwirepl(j ,iw(j)) 

enddo 
enddo 
enddo 
enddo 
enddo 
enddo 
enddo 
enddo 
enddo 
endif 

if (ncomb.gt . 18000) then ! if too many combinations, skip it 

write (6,*) 'ncomb=' ,ncomb 

return 
endif 



C get chisq and pointing in v for each combination 
chimin=99999. 
pointmin=9999 . 
DO M=l,ncomb 
Szz=0 
Sz=0 
Sv=0 
Szv=0 

DO j=l,nchain 

k=iichcmi( j ) 

u(k)=a_u(i)+b_u(i)*Z(k) 

v(k)=(u(k) *sth(k) -UO (k) *sth(k)+VO (k) *cth(k)+ 

+ wire(M, j)*wiresp)/cth(k) 
Szz=Szz+Z(k)**2 
Sz=Sz+Z(k) 
Sv=Sv+v(k) 
Szv=Szv+Z(k)*v(k) 
ENDDO 

denom=nchcim*Szz-Sz**2 
av= (Szz*Sv-Sz*Szv) / denom 
bv= (ncham*Szv-Sz*Sv) / denom 
chi(M)=0 
DO J=l,ncham 
K=iicham( J) 
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predwire= (- (u (k) -UO (K) ) *siii (ang (K) ) 
+ + (av+bv*Z (K) -VO (K) ) *cth (K) ) /wiresp 

chi (M) =chi (M) + (predwire-wire (M , J) ) **2 
ENDDO 

chi (M) =chi (M) / (sigsq* (ncham-2) ) 
if (chi (M) . It . chimin) then 

chiinin=chi (M) 

ichi=M 
endif 

thtail=ataii(bv) 
vlead=av+bv*149.8 

thgain=atan(vlead/ (149 . 8+a_u(i) /b_u(i) ) ) 
sigth=sqrt ( (0 . 01*sigsq) *ncham/ (ncham*Szz-Sz**2) ) * 
+ cos(thtail)**2 

point (M) =abs (thgam-thtail) /sigth 
if (point (M) . It .pointmin) then 

pointmin=point (M) 

ipoint=M 
endif 
ENDDO 

C require vpoint < 8 and chisq < 7, give preference to best pointing 
if (ncomb.gt.O) then 
icomb=0 

if ((pointmin.lt.8.) -and. (chi(ipoint) .lt.7)) then 

icomb=ipoint 
else 

if ((chimin.lt.7) .and. (point(ichi) .lt.8)) then 
icomb=ichi 
endif 
endif 

nwire=nnwire (i) 
nnwire(i)=0 

if (icomb.gt.O) then !no good track — > nnwire(i)=0 

ichlast=0 

DO J=l,nwire ! refill list of wires corresp 

ichaiii=int((wires(I, J)-l)/128.)+l ! to track 
if (uchamdcham) . eq. 1) then 

nnwire(i)=nnwire(i)+l 

wires (i ,nnwire (i) ) =wires (I , J) 
else 

if (ichcun.ne . ichlast) then 
do k=l,ncham 
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if (iichamCk) . eq. icham) then 
nnwire ( i ) =iiiiwire ( i ) +1 

wires (i .nnwire (i) ) = (icham-1) *128+wire (icomb , k) 
endif 
enddo 
endif 
endif 

ichlast=icham 
ENDDO 
endif 
endif 



999 continue 
enddo 

C refill wires. inc 
ntr=ntrack 
ntrack=0 
nchl=0 
do i=l,ntr 

if (nnwire (i) .gt.O) then 
ntrack=ntrack+l 
if (i.le.nch) nchl=nchl+l 
nnwire (ntrack) =nnwire ( i) 
do j=l, nnwire (i) 

wires (ntrack , j ) =wires (i , j ) 
enddo 
endif 
enddo 

C refit tracks from new wires. inc (for all tracks is cpu intensive but easier) 
call fit_wires 



RETURN 
END 



Appendix C 
Vertexer code 



C.l Vertexer: vertexer. f 

SUBROUTINE vertexer 
IMPLICIT NONE 

INTEGER HLIM.HBSIZE 

PARAMETER(HBSIZE=500000) 
COMMON/PAWC/HLIM (HBSIZE) 

integer i, j ,k,l,in,keep 

real Ml (50, 2, 2) ,M1INV(50,2,2) ,MM(2,2) ,MMINV(2,2) 
real dchisq,iimean,vinean,denoin 

integer nmatch, iinatch(50) , isinatched(lOO) 

integer nnclust(50) ,iiclust(50,50) , iclust ,new,newtrack,ninax,nncl 
integer dstnumb(lOO) ,nv 

integer max, imax,nsame ,nscimeu,maxu, icham, ichlast 

integer nback(lOO) ,nb(24) 

real thtail , sigth , thgam , zlead , siguO 

integer ntr ,nhead,ntail 

real upoint , vpoint , upoint2 , vpoint2 

integer nch.l,nch.2 
common/newnch/nchl ,nch2 
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real uoldl ,voldl 

real dumin , dvmiii , dmin , du , dv , uold , void 

common/t aildist/ dumin , dvmin , dmin , du , dv , uold , void 

include 'fit.inc' 
include ' wires. inc' 
include 'vrtx.inc' 
include 'VRTCS.INC 

include ' track_dat . inc ' 
include 'dst.inc' 

zlead=149 . 82 

C fill wires. inc with all track segments 
ntrack=0 

DO k=l,nch+nnu 

C get number of wires (don't distinguish betwen heads or tails) 
IF (k.le.nch) THEN 

ntrack=ntrack+l 
if (ntrack.gt.lOO) write(6,*)ntrack 
dstnumb (ntrack) =k 
nnwire (k)=int (charged (k, 6) ) 
DO J=l , nnwire (ntrack) 

wires (k , j ) =wires_c (k , J) 
ENDDO 
ELSE 

DO L=9,24 
nb(L)=0 
ENDDO 

DO j=l, neutral (k-nch, 6) 

ichcim=int ( (wires_n(k-nch, j)-l)/128. )+l 

nb(icham)=l 
ENDDO 

nback (ntrack+1) =0 
DO L=9,24 

nback (ntrack+1 ) =nback (ntrack+1 ) +nb (L) 

ENDDO 

if (nback (ntrack+1) .ge. 8) then lean chainge to exclude mid-8 tracks 
ntrack=ntrack+l 

dstnumb (ntrack) =k 

nnwire (ntrack) =int (neutral (k-nch , 6 ) ) 
DO J=l , nnwire (ntrack) 
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wires (ntrack , j ) =wires_n (k-nch , J) 

ENDDO 

endif 
ENDIF 
ENDDO 

C fit tracks 

call fit_wires 
call refit 

C get position of tracks at lead and uncertainty in pos 
call vrtx(zlead) 

C remove heads which don't point and tracks which don't hit lead 
ntr=ntrack 
ntrack=0 
nch2=0 
do i=l,ntr 

if (i.le.nchl) then 

upoint=-a_u ( i ) /b_u ( i ) 
else 

upoint=0 
endif 

C for now, don't demand tails point, loose on heads 

if ( (upoint .gt . -50) . Eind. (upoint . It . 60 . ) ) then 
if (i.le.nchl) then 
thtail=atan(b_v(i) ) 

thgain=atan(v(i)/ (zlead+a_u(i)/b_u(i) ) ) 

sigth=sqrt (covar (i ,4,4)) *cos (thtail) **2 

vpoint=abs (thgam-thtail) / sigth 
else 

vpoint=0 
endif 

if ( (vpoint . It . 8 . ) .and. (chisq(i) . It. 20. ) .and. 
+ (u(i) .gt.3.3) .and. (u(i) .It. 11. 3) .and. 

+ (v(i) .gt.-4.0) .and. (v(i) .lt.4.0)) then 

C refill wires. inc with kept track segments 
ntrack=ntrack+l 
dstnumb (ntrack) =dstnumb(i) 

if (i.le.nchl) nch2=nch2+l ! number of heads which point 
nnwire (ntrack) =nnwire ( i) 
do j=l , nnwire (i) 

wires (ntrack , j ) =wires (i , j ) 
enddo 
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a_u (nt rack) =a_u ( i ) 
b_u(ntrack)=b_u(i) 
a_v(ntrack)=a_v(i) 
b_v (ntrack) =b_v (i) 
do j=l,4 
do k=l,4 

covar (ntrack, j ,k)=covar (i , j ,k) 
enddo 

enddo 

chisq(ntrack)=chisq(i) 
u (ntrack) =u(i) 
v(ntrack)=v(i) 
siguu (ntrack) =siguu(i) 
siguv (ntrack) =siguv(i) 
sigvv (ntrack) =sigvv(i) 

endif 

endif 
enddo 



C 

C get i,j for tracks to match 
nvrtcs=0 



C increase uncer in pos of tails to account for multiple scattering 
do i=l, ntrack 
ismatched(i)=0 
if ( i . gt . nch2 ) then 

if (siguu(i) .It. 0.007) siguu(i)=.007 
if (sigvv(i) .It. 0.092) sigvv(i)=.092 
endif 
enddo 



C first match heads to other segments (smaller uncertainty) 
do i=l, ntrack 

if (i.gt.nch2) goto 555 
nmatch=0 
do j=l, ntrack 

if (j.eq.i) goto 557 

CALL match(i, j ,dchisq) 

if ((dchisq/2.) .lt.5.) then 

nmatch=nmatch+l ! number of matches to segment i 

imat ch (nmat ch) = j 

ismatched(j)=l 
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endif 
557 continue 

enddo ! do j 

iniatch(nmatch+l)=i ! (include track i in list of matched tracks) 



C if this head has already been found to match other segments, 
C combine all matches 
new=l 

do k=l,nvrtcs 

do l=l,nnclust(k) 

do j=l,nmatch+l 

if (imatch(j) .eq.iiclust(k,l)) then 
new=0 
iclust=k 

endif 
enddo 
enddo 
enddo 

if (new.eq.l) then 

nvrt c s=nvrt c s+ 1 
iclust=nvrtcs 
nnclust (iclust)=0 
endif 

nmax=nnclust (iclust) 
do j=l,nmatch+l 
newtrack=l 
do l=l,nmax 

if (imatch(j) .eq.iiclust(iclust,l)) newtrack=0 
enddo 

if (newtrack. eq. 1) then 

nnclust (iclust)=nnclust(iclust)+l 
iiclust (iclust , nnclust (iclust) )=imatch(j ) 

endif 

enddo 

555 continue 

enddo !do i 



C now match tail-tail to get gammas 
uoldl=999. 
do i=l,ntrack 

if ((i.le.nch2) .or. (ismatched(i) .eq.l)) goto 556 
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nmatch=0 
uoldl=u(i) 
voldl=v(i) 
dmiii=999 . 
duiiiin=999 . 
dvmin=999 . 
do j=l,ntrack 

if ((j .gt.iich2) .and. (j .ne.i) .and. (ismatchedCj) .ne.l)) then 

CALL match(i , j , dchisq) 

if ((dchisq/2.) .lt.5.) then 

nmat ch=nmat ch+ 1 

imat ch (nmat ch) = j 
endif 

endif ! j ne i 
enddo ! do j 

imat ch (nmat ch+ 1 ) =i 

new=l 

do k=l,nvrtcs 

do l=l,nnclust(k) 

do j=l,nmatch+l 

if (imatch(j) .eq.iiclust(k,l)) then 
new=0 
iclust=k 

endif 
enddo 
enddo 
enddo 

if (new.eq.l) then 

nvrtcs=nvrtcs+l 

iclust=nvrtcs 

nnclust (iclust)=0 
endif 

nmax=nnclust (iclust) 
do j=l,nmatch+l 
newtrack=l 
do l=l,nmax 

if (imatch(j) .eq.iiclust (iclust,!)) newtrack=0 
enddo 

if (newtrack. eq. 1) then 

nnclust (iclust)=nnclust(iclust)+l 

iiclust (iclust , nnclust (iclust) )=imatch(j ) 
endif 
enddo 
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continue 
enddo !do i 

if (uoldl.lt. 11.25) then 

uold=uoldl 
vold=voldl 
endif 



fill VRTCS.INC with found vertices 
nv=0 

do k=l,nvrtcs 
nverts (k, 1)=0 
nverts(k,2)=0 
do j=l .nnclust (k) 
i=iiclust(k, j) 
if (i.gt.nch2) then 

nverts (k , 2) =nverts (k , 2) +1 
else 

nverts (k, l)=nverts(k, 1)+1 
endif 
enddo 

if ((nverts(k,l) .eq. 1) .and. (nverts(k,2) .gt.O)) goto 123 

! if (i,n>0) charged track, pointing is already ok 

nncl=nnclust (k) 
nnclust (k)=0 
nhead=nverts (k, 1) 
ntail=nverts (k , 2) 
nverts (k,2)=0 
nverts (k,l)=0 
do j=l,nncl 

i=iiclust(k, j) 

upoint=-a_u(i) /b_u(i) 

thtail=atan(b_v(i)) 

thgam=atan(v(i) / (zlead+a_u(i) /b_u(i) ) ) 

sigth=sqrt(covar(i,4,4))*cos(thtail)**2 

vpoint=abs (thgam-thtail) /sigth 

if ( ( (upoint . gt . -40) . Eind. (upoint .It . 50) . and. 
+ (vpoint.lt. 8)) .or. 

+ ((nhead.eq.O) .and. (ntail.gt.l) .and. 
+ (upoint. It. 50))) then 
keep=l 

if ((nhead.eq.O) .and. (ntail.eq.l)) then Itighter on (0,1) 's 
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keep=0 
do 1=9,24 
nb(l)=0 
enddo 

do 1=1 ,niiwire(i) 

icham=int ( (wires (i , 1) -1)/128 . )+l 

nb(ichain)=l 
enddo 
ichcmi=0 
do 1=9,24 

i cham= i cham+nb ( 1 ) 
enddo 

if (icham.ge. 14) keep=l 
endif 

if (keep.eq.l) then 
nnclust (k) =nnclust (k) +1 
iiclust (k , nnclust (k) )=i 
if (i.gt.nch2) then 

nvert s (k , 2) =nverts (k , 2) +1 
else 

nvert s (k, l)=nverts (k, 1)+1 
endif 
endif 
endif 
enddo 
123 continue 

if (nnclust (k) .gt . 0) then 
nv=nv+l 

nverts(nv, l)=nverts(k, 1) 
nvert s (nv , 2) =nverts (k, 2) 
vrtxpos (nv , 1) =0 
vrtxpos (nv , 2) =0 
if (nnclust (k) .gt. 20) write (6, *) 'nnclust=' , nnclust (k) 
C get mean position of vertex 
MM(1,1)=0 
MM(1,2)=0 
MM(2,2)=0 

do j=l , nnclust (k) 
i=iiclust(k, j) 
nvert s (nv , j +2 ) =dstnumb ( i ) 
nverts (nv , j +52) =i 
Ml(j,l,l)=siguu(i) 
Ml(j,2,2)=sigvv(i) 
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Ml(j,l,2)=siguv(i) 
Ml(j,2,l)=Ml(j,l,2) 

denom=Ml(j ,l,l)*Ml(j ,2,2)-Ml(j ,l,2)*Ml(j ,2,1) 
MlINV(j,l,l)=Ml(j,2,2)/denoin 
MlINV(j,2,2)=Ml(j ,l,l)/denoin 
MlINVCj ,1,2)=-Ml(j ,l,2)/denom 
MIINV ( j , 2 , 1) =-Ml ( j , 2 , 1) /denom 
MM(l,l)=MM(l,l)+MlIMV(j,l,l) 
MM(l,2)=MM(l,2)+MlINV(j ,1,2) 
MM(2,2)=MM(2,2)+MlINV(j,2,2) 
enddo 

MM(2,1)=MM(1,2) 

denom=MM(l,l)*MM(2,2)-MM(l,2)*MM(2,l) 

MMINV ( 1 , 1 ) =MM ( 2 , 2 ) /denom 

MMINV ( 2 , 2 ) =MM ( 1 , 1 ) /denom 

MMINV ( 1 , 2 ) =-MM (1,2) /denom 

MMINV ( 2 , 1 ) =-MM ( 2 , 1 ) /denom 

umean=0 

vmean=0 

do j=l,nnclust(k) 
i=iiclust(k, j) 

umean=umean+MMINV (l,l)*MlINV(j,l,l)*u(i)+ 



+ MMINV(l,2)*MlINV(j ,2,l)*u(i)+ 

+ MMINV(l,l)*MlINV(j ,l,2)*v(i)+ 

+ MMINV(l,2)*MlINV(j ,2,2)*v(i) 

vmean=vmean+MMINV (2 , 1 ) *M1 INV ( j , 1 , 1 ) *u ( i ) + 

+ MMINV(2,2)*MlINV(j ,2, l)*u(i)+ 

+ MMINV(2,l)*MlINV(j ,l,2)*v(i)+ 

+ MMINV(2,2)*MlINV(j ,2,2)*v(i) 



enddo 

vrtxpos(nv, l)=umeEin 
vrtxpos (nv , 2) =vmean 
endif 

enddo 

nvrtcs=nv 

RETURN 
END 

C 

C get u,v,uncer at lead 

SUBROUTINE vrtx(zlead) 
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include 'fit.inc' 
include ' wires. inc' 
inc lude ' vrtx . inc ' 

include ' track_dat . inc ' 
include 'dst.inc' 

thu=45 . 43*3 . 141592654/180 . 

cu=cos (thu) 

su=sin(thu) 

do i=l,ntrack 
u(i)=a_u(i)+b_u(i)*zlead 

v(i)=a_v(i)+b_v(i)*zlead 

siguu(i)=covar(i,l,l)+2*zlead*covar(i,l,2)+zlead**2*covar(i,2,2) 
sigvv (i) =covar (i ,3,3) +2*zlead*covar ( i , 3 , 4) +zlead**2*covar (i , 4,4) 
siguv ( i ) =covar ( i , 1 , 3) +zlead* (covar (i ,1,4) +covar (i ,2,3))+ 
+ zlead**2*covar(i,2,4) 
enddo 

RETURN 
END 

C 

C get chisq of match between pair of tracks 

SUBROUTINE matchCi, j ,dchisq) 

IMPLICIT NONE 

include 'fit.inc' 
include ' vrtx . inc ' 

integer i,j 

real CI (2,2) ,C2(2,2) ,C1INV(2,2) ,C2INV(2,2) ,C(2,2) ,CINV(2,2) 
real dchisq,uinean,vineaii,denoin 

integer nchl,nch2 
common/newnch/nchl ,nch2 

real dist 

real dumin , dvmin , dmin , du , dv , uold , void 
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common/t aildist/ dumin , dvmin , dmin , du , dv , uold , void 

real upoint , vpoint ,upoint2 , vpoint2 
real thtail , sigth,thgam,zlead,siguO 

INTEGER HLIM.HBSIZE 
PARAMETER (HBSIZE=500000) 
COMMON/PAWC/HLIM (HBSIZE) 

Cl(l,l)=siguu(i) 

Cl(2,2)=sigvv(i) 

Cl(l,2)=siguv(i) 

C1(2,1)=C1(1,2) 

C2(l,l)=siguu(j) 

C2(2,2)=sigvv(j) 

C2(l,2)=siguv(j) 

C2(2,1)=C2(1,2) 

C invert 

denom=Cl(l,l)*Cl(2,2)-Cl(l,2)*Cl(2,l) 

ClINV(l,l)=Cl(2,2)/denoin 

ClINV(2,2)=Cl(l,l)/denom 

ClINV(l,2)=-Cl(l,2)/denom 

C1INV(2 , 1) =-Cl (2 , 1) /denom 

denom=C2 ( 1 , 1) *C2 (2 , 2) -C2 (1 , 2) *C2 (2 , 1 ) 

C2INV ( 1 , 1) =C2 (2 , 2) /denom 

C2INV(2,2)=C2(l,l)/denoin 

C2INV ( 1 , 2) =-C2 (1,2) /denom 

C21NV (2 , 1) =-C2 (2 , 1) /denom 

C(1,1)=C1INV(1,1)+C2IMV(1,1) 

C(2,2)=C1INV(2,2)+C2INV(2,2) 

C(1,2)=C11NV(1,2)+C21NV(1,2) 

C(2,1)=C(1,2) 

denom=C(l,l)*C(2,2)-C(l,2)*C(2,l) 
CINV ( 1 , 1 ) =C (2 , 2) /denom 
CINV (2 , 2) =C (1 , 1) /denom 
CINV ( 1, 2) =-C( 1,2) /denom 
CINV (2, 1)=-C (2,1) /denom 

umean=CINV(l,l)*(ClIMV(l,l)*u(i)+C2IMV(l,l)*u(j))+ 



+ CIMV(l,2)*(ClINV(2,l)*u(i)+C2INV(2,l)*u(j))+ 
+ CINV(l,l)*(ClINV(l,2)*v(i)+C2INV(l,2)*v(j))+ 
+ CINV(l,2)*(ClIMV(2,2)*v(i)+C2INV(2,2)*v(j)) 
vmean=CINV(2,l)*(ClINV(l,l)*u(i)+C2INV(l,l)*u(j))+ 
+ CINV(2,2)*(ClINV(2,l)*u(i)+C2INV(2,l)*u(j))+ 
+ CINV(2,l)*(ClIMV(l,2)*v(i)+C2INV(l,2)*v(j))+ 
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+ CINV(2,2)*(ClINV(2,2)*v(i)+C2INV(2,2)*v(j)) 

dchi sq= (u ( i ) -umecin) *C 1 INV (1,1)* (u(i) -umecin) 

+ +2*(u(i)-umean)*ClINV(l ,2)*(v(i)-vmean) 

+ +(v(i)-vmeaii)*ClINV(2,2)*(v(i)-vmean) 

+ +(u(j)-umean)*C2INV(l,l)*(u(j)-umean) 

+ +2* (u ( j ) -umean) *C2INV (1 , 2) * (v ( j ) -vmean) 

+ + (v ( j ) -vmean) *C2INV (2 , 2) * (v( j ) -vmean) 



RETURN 
END 



C 

C vrtx . inc 

real u(200) , v(200) , siguu(200) ,sigvv(200) ,siguv(200) 
common/ lead/u , v , siguu , sigvv , siguv 

C (determined by subroutine vrtx in vertexer.f) 
C position eind uncertainty of track at z=zlead 

C 

C VRTCS . INC 

integer nvrtcs 

integer nverts(100, 102) 

real vrtxpos(100,2) 

integer ischarged(lOO) , isgamma(lOO) 

common/fndvrts/nvrtcs ,nverts, vrtxpos , ischarged, isgaimna 

C (determined by vertexer.f) 

C nvrtcs=number of vertices at zlead 

C nverts (i , l)=number of heads in vertex i 

C nverts(i,2)=number of tails in vertex i 

C nverts (i,2+j)=dst-track-number of track segment j in vertex i 

C nverts(i,52+j)=f itter-track-number of track segment j in vertex i 



C (determined by chgam.f) 

C ischarged (i)/isgamma(i) =1 if vertex is a charged track/ gamma 
C =0 otherwise 
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C.2 Code for classifying vertices as charged 
tracks or photon conversions: chgam.f 

SUBROUTINE CHGAM 

IMPLICIT NONE 

include 'EVENT. INC 

include 'fit.inc' 

include 'wires. inc' 

include 'VRTCS.INC 

integer nch.l,nch2 
common/ newnch/ nchl , nch2 
INTEGER NN 

REAL sigsq,wiresp,pi ,ucing, sum, cu, su 
REAL U0(24) ,V0(24) ,ang(24) ,Z(24) 

common/ cham_parain/NN ,sigsq,wiresp,pi, uang , sum , cu , su , UO , VO , ang , Z 

INTEGER I, J, IDS, k, ihit , icham, nb(24), keep.ii, j j ,kk,nu,nv 
INTEGER nkeep,ikeep(20) ,imin 
REAL dist,dmin,r 

real mean, sigmean,upoint , sigupoint ,wire 

integer itr , itrO ,keepcham,keepwire ,keeppoint ,npoint 



do j=l,nvrtcs 
ischarged(j)=G 
isgamma(j)=0 



C Cyrus's acceptance 

if ((vrtxposCj ,1) .gt.4.25) .and. (vrtxposCj ,1) .It. 10.25) .and. 
+ (vrtxpos ( j ,2) .gt . -3 . ) . and . (vrtxposC j , 2) . It .3. ) ) then 

r=sqrt( (vrtxposCj , l)-7.25)**2+vrtxpos(j ,2)**2) 
if (r.lt.4.) then 

C (l,n>0) charged track 

if ((nvertsCj , 1) .ge . 1) .and. (nvertsCj ,2) .gt .0)) ischargedCj )=1 

C (1,0) 's have to pass tighter cuts 

if ( (nverts(j , 1) .ge. 1) .and. (nverts(j ,2) .eq.O)) then 
keep=0 



C require hit all 8 chambers 
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do i=l,nverts(j ,1) 
itrO=nverts(j ,52+i) 

keepch£mi=0 

do k=l,8 
nb(k)=0 

enddo 

do k=l,nnwire(nverts(j ,52+i)) 
nbCintC (wires (nverts(j,52+i) ,k)-l)/128. )+l)=l 
enddo 
icheim=0 
do k=l,8 

if (nb(k) .eq.O) then 

wire=nint ( (- (a_u(itrO) +b_u(itrO) *Z(K) -UO (K) ) *sin(aiig(K) ) 
+ +(a_v(itrO)+b_v(itrO)*Z(K)-VO(K))*cos(ang(K)))/wiresp ) 

dmin=128 

do jj=l,nhits 

if ((inod(chamb_nijm(jj) ,100) .eq.k) .and. 
+ (abs(wire_num(j j)-wire) .It.dmin)) 

+ dinin=abs (wire_nuin(j j )-wire) 

enddo 

if (dmin.le.l) then 
nb(k)=l 

nnwire (itrO)=nnwire(itrO)+l 
wires (itrO , nnwire (itrO) ) =wire+ (k-1) *128 
endif 
endif 
ichcmi=icham+nb (k) 
enddo 

if (icham. eq. 8)keepcham=l 

C don't share wires in all 3 u or all 3 v cheimbers with other head 
keepwire=l 
do jj=l,nvrtcs 
if (jj .ne. j) then 
do ii=l ,nverts(j j , l)+nverts(j j ,2) 
if (nvertsCj j ,52+ii) .le.nch2) then 
itr=nverts ( j j ,52+ii) 
nv=0 
nu=0 

do kk=l , nnwire (itr) 
do k=l , nnwire (itrO) 

if (wires(itr,kk) .eq.wires(itrO,k)) then 
ichain=int ( (wires (itrO , k) -1) /128 . ) +1 
if ((icham. eq.l) .or. (icham. eq. 3) .or. (icham. eq. 6) ) 
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+ nv=nv+l 

if ((icham.eq.2) .or. (icham.eq.5) .or. (icham.eq.8)) 
+ nu=nu+l 
endif 
enddo 
enddo 

if ( (nv . eq . 3) . or . (nu . eq . 3) ) keepwire=0 

endif 

enddo ! ii 

endif 

enddo ! j j 

C point to mean z of other charged tracks 
if (i.eq.l) then 
mean=0 
sigmeaii=0 
npoint=0 
do jj=l,nvrtcs 

if ( (nverts (j j , 1) . ge . 1) . and. (nvertsCj j ,2) . ge . 1) ) then 
do ii=l , nverts (jj , l)+nverts(j j ,2) 
if (nvertsCj j ,52+ii) .le.nch2) then 

npoint=npoint+l 

itr=nverts(j j ,52+ii) 

upoint=-a_u(itr) /b_u(itr) 

sigupoint= ( covar (itr ,1,1)+ (a_u (itr) /b_u (itr) ) **2* 
+ covar(itr,2,2)-2*(a_u(itr)/b_u(itr))*covar(itr,l,2) ) 
+ /b_u(itr)**2 

mean=mean+upoint/sigupoint 
sigmean=siginean+l . / sigupoint 
endif 
enddo ! ii 
endif 
enddo ! j j 

if (npoint .gt .0) then 

sigmean=l . /sigmean 

me an=me cin* s i gme an 
endif 

endif ! i=l 

keeppoint=l 

if (npoint .gt .0) then 

dist=((-a_u(itrO)/b_u(itrO))-mean) 
+ /sqrt (sigmean+( covar (itrO , 1 , l)+(a_u(itrO) / 

+ b_u (itrO) ) **2*covar (itrO , 2 , 2) -2* (a_u (itrO) / 

+ b_u(itr0))*covar(itr0,l,2) )/b_u(itr0)**2 ) 
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if (abs (dist) . gt . 2 . ) keeppoint=0 
endif !npoint>0 

if ( (keepcham. eq. 1) . Euid. (keepwire . eq. 1) . eind. 
+ (keeppoint . eq. 1) ) keep=l 

enddo ! i 

if (keep.eq.l) ischarged(j)=l 
endif ! (n,0) 
C (0,n) gamma 

if ((nvertsCj ,1) .eq.O) .and. (nvertsCj ,2) .gt.O)) isgammaCj )=1 

endif ! vrtxpos 
endif !r<4 
enddo ! j 



RETURN 
END 
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